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ABSTRACT 



Context. The identification and study of the first galaxies remains one of the most exciting topics in observational cosmology. The 
determination of the best possible observing strategies is a very important choice in order to build up a representative sample of 
spectroscopically confirmed sources at high-z (z > 7), beyond the limits of present-day observations. 

Aims. This paper is intended to precisely adress the relative efficiency of lensing and blank fields in the identification and study of 
galaxies at 6; z; 12. 

Methods. The detection efficiency and field-to-field variance are estimated from direct simulations of both blank and lensing fields 
observations. Present known luminosity functions in the UV are used to determine the expected distribution and properties of distant 
samples at z i 6 for a variety of survey configurations. Different models for well known lensing clusters are used to simulate in details 
the magnification and dilution effects on the backgound distant population of galaxies. 

Results. The presence of a strong-lensing cluster along the line of sight has a dramatic effect on the number of observed sources, 
with a positive magnification bias in typical ground-based "shallow" surveys (AS; 25.5). The positive magnification bias increases 
with the redshift of sources and decreases with both depth of the survey and the size of the surveyed area. The maximum efficiency 
is reached for lensing clusters at z ~ 0.1 - 0.3. Observing blank fields in shallow surveys is particularly inefficient as compared to 
lensing fields if the UV LF for LBGs is strongly evolving at z i 7. Also in this case, the number of z > 8 sources expected at the typical 
depth of JWST (AB ~ 28 - 29) is much higher in lensing than in blank fields (e.g. a factor of ~ 10 for AB < 28). All these results have 
been obtained assuming that number counts derived in clusters are not dominated by sources below the limiting surface brightness of 
observations, which in turn depends on the reliability of the usual scalings applied to the size of high-z sources. 
Conclusions. Blank field surveys with a large field of view are needed to prove the bright end of the LF at z i 6 — 7, whereas lensing 
clusters are particularly useful for exploring the mid to faint end of the LF. 

Key words, galaxies : formation - galaxies : high redshift - galaxies : photometry - galaxies : clusters : lensing - 



■ »-h , 1 . Introduction 



Constraining the abundance of z > 7 sources remains an 
important challenge of modern cosmology. Recent WMAP re- 
sults place the first building blocks of galaxies at redshifts 
z = 11.0 + 1.4, suggesting that reionization was an ex- 
et al.ll2009l) . 



tended process (Dunkley < 



). Distant star-forming sys- 



tems could have been responsible for a significant part of the 
cosmic reionization. Considerable advances have been made 
during the last years in the observation of the early Universe 
with the discovery of galaxie s at ~6-7, close to the end 
of the reionization epoch (e.g. [Hu et al] )2002t iKodaira et al] 
20031; ICubv et al.ll2003t Ikneib et al.ll20'04t |stanwav et al.[|2004|: 



Bouwens et al. 2004al 120061: live et al.H2006UBradlev e t al. 2008; 
Zhen g et al.l2009l), and the first prospects up to z~10 (|Pello et al] 
2004t IRichard et alj|2006tlStark et alj2007t iRichard et alT2 008 ; 
Bouwens et al.ll2008ll2009d) . 



High-z surveys are mainly based on two different and com- 
plementary techniques: the dropout (Lyman-o' Break) identifica- 
tion, which is an extrapolation of the drop-out technique used 
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for Lyman-Br eak Galaxies (LBGs.lSteidel et alJI 19991) to highe r 
redshifts (e.g. iBouwens et alj|2006h IRichard et al.ll2006l |2008), 
and the narrow-band (N B) imaging aim ed at detecting Lyman 
a emitters (LAEs, e.g . [Taniguchi et al] 120051: live et alJ 12 006: 
iKashikawa etaT] 12000: Willis et al] 120061: ICubv et alj 120071). 
Using the former technique. IBouwens et al.l d2008l) found a 
strong evolution in the abundance of galaxies between z~7-8 and 
z~3-4, the SFR density beeing much smaller at very high-z up 
to the limits of their survey, in particular towards the bright end 
of the Luminosity Function (LF). A strong evolution is also ob- 
served in the number density of sources detected with NB tech- 
niques, which see ms to be much smaller at z > 7 than in the 
z~ 5 - 7 interval dive et alj|20"06l ICubv et al.|[2007t IWillis et al.l 
120081) . 

Both dropout and NB approaches require a subsequent spec- 
troscopic confirmation of the selected candidates. For now ap- 
proximately ten galaxies b eyond z~6.5 are known with se- 
cure spectroscopic redshifts (Hu et al. 2002; Kodaira et al. 2003; 
ICubv et all 120031: iTaniguchi etaT] 120051: llveetal.1 120061) . All 
samples beyond this redshift are mainly supported by photo- 
metric considerations dKneib et al]l2004t IBouwens et al] l2004b. 
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2006- lRichard"eFaIll2006l 120081: iBradlev et al.ll2008l) . This situ- 
ation is expected to dramatically improve in the near future with 
the arrival of a new generation of multi-object spectrographs 
in the near-IR, such as MOIRCS/Subaru, Flamingos2/Gemini- 
S (~ 2009), or EMIR/GTCQ (~ 2012), with well suited field 
of view, spectral resolution and sensitivity. These forthcoming 
facilities should provide spectroscopic confirmation for a large 
number of z>7 candidates identified from deep photometric sur- 
veys, and the first characterization of the physical properties of 
these sources (e.g. IMF, stellar populations, fraction of AGN, ...). 

The aim of this paper is to determine the best possible ob- 
serving strategies in order to build up a representative sam- 
ple of spectroscopically confirmed z > 7 galaxies. The photo- 
metric pre-selection of high-z candidates could be achieved ei- 
ther in blank fields or in lensing clusters. This later technique, 
also first referred to as the "gravitational telescope" by Zwicky 
, has proven highly successful in identifying a large fraction 
of the most distant galaxies known today thanks to magnifi- 
cations by typically 1-3 magnitudes in the cluster core (e.g. 
Ellis et al.ll200lHHu et al.l2002l;lKneib et alJ2004tlBradlev et alj 
20081: IZheng et alj|2009t Bradac et al.1120091) . The presence of a 
strong lensing cluster in the surveyed field introduces two op- 
posite effects on number counts as compared to blank fields. In 
one hand, gravitational magnification increases the number of 
faint sources by improving the detection towards the faint end 
of the LF. On the other hand, the reduction of the effective sur- 
face by the same magnification factor leads to a dilution in ob- 
served counts. The global positive/negative magnification bias 
obvio usly depends on the slop e of the number counts, as pointed 
out bv lBroadhurst et all d 1995b . 

This paper addresses the relative efficiency of surveys con- 
ducted on blank and lensing fields as a function of the relevant 
parameters, namely the redshift of the sources, the redshift and 
properties of the lensing clusters and the survey characteristics 
(i.e. area, photometric depth...). This calculation requires a de- 
tailed simulation of observations using lensing models, and re- 
alistic assumptions for the properties of background sources ac- 
cording to present-day observational results, in particular for the 
luminosity function and typical sizes of z > 7 galaxies. 

The paper is organized as follows. In Section|2]we describe 
the simulations performed in order to determine the relative de- 
tection efficiency for high-z sources, both in lensing and blank 
fields. Section [3] presents the results, in particular the detection 
efficiency achieved as a function of redshift for both sources and 
lensing clusters, together with a discussion on the influence of 
lensing cluster properties and field-to-field variance. A discus- 
sion is presented in Section|4]on the relative efficiency as a func- 
tion of survey parameters, and a comparison between simula- 
tions and present surveys. Conclusions are given in Section[5] 

Throughout this paper, we adopt a concordance cosmo- 
logical model, with £2a = 0.7, Q„, = 0.3, and Hq 
70 km s Mpc~ l . All magnitudes are given in the AB system. 
Conversion values between Vega and AB systems for the filters 
used in this paper are typically Cab- 0.95, 1.41 and 1.87 in J, H 
and K respectively, with m AB = m Vega + Cab- 



2. Simulations of lensing and blank field 
observations 

2.1. Simulation parameters 

This section describes the ingredients used in the simulations to 
implement different assumptions that would affect our efficiency 
in detecting high redshift galaxies. There are three important as- 
pects to be considered in the comparison between lensing and 
blank fields. The first one is the LF and typical sizes of sources. 
The second one concerns the properties of the lensing clusters, 
in particular their mass distribution and redshift. The third one is 
related to the survey parameters, namely the photometric depth 
and the size of the field. All these aspects are discussed in this 
section. Table[TJprovides the list of parameters used in these sim- 
ulations, together with the range of values explored. 

2.1.1. Source Properties 

These simulations are focused on the detection of sources in the 
redshift range 6 < z < 12, a relevant domain for spectroscopic 
follow-up with near-infrared instruments. The lower limit of this 
redshift domain overlaps with c urrent photometric su rveys mea- 
suring the LF at z ~ 6 - 7 (e.g. lBouwens et al. 2007). However, 
the LF is still largely unconstrained beyond z Z 7 because of the 
lack of spectroscopic confirmation of photometric samples and 
the relatively small size of the surveyed volumes. 

The abundance of background sources at these redshifts is 
given by the luminosity function <p(L), with L the rest-frame 
UV luminosity at 1500 A. cp(L) is the most basic description 
of the galaxy population from an observer point of view. We 
adopt a parametriza tion based on the analytical Schechter func- 
tion (ISchechterll 19761) : 

cf>(L)dL = f^ expf--^)^) (1) 

The slope at faint luminosities a, the characteristic luminos- 
ity L* and the normalization factor O* have been constrained 
by several photometric survey s targeting Lyman Break Galaxies 
(LBG) at high redshift (z > 4) dOuchi et al.l2004irBeckwith et al 
| 2006l IBouwens etal. 1 120071 iBouwens et alj|2008t iMcLure et al. 
2009, iHenrv et all 120091) . Three different representative cases 
are being discussed in our simulations, basically exploring our 
present knowledge (or lack of knowledge) of the overall shape 
oftheLFatz>6(Table[B: 

(a) An "optimistic" scenario where LBGs show no-evolution 
from z ~ 6, with th e LF parameters as determined by 
iBeckwith etail (120061) . Indeed, the LF at z ~ 6 found 
by these authors d isplay the same shape as for z ~ 3 
( Stei del et al.ll999|), but a 3 times sm aller normalization fac- 
tor (but see, e.g. IMcLure et al.ll2009l) . 

(b) A constant L F based on robust measurements by 

IBouwens et ail (120071) at z ~ 6 in the Hubble UDF 
but usi ng the more recent fit parameters from lBouwens et ail 
(2008). As compared to model (a), this LF exhibits a 
turnover towards the bright end. 

(c) The ev olutionary LF recently proposed by IBouwens et alj 
d2008), which includes an important dimming of L* with in- 
creasing redshift. This LF represents the "pessimistic" case 
with respect to the case (a), with very few high-luminosity 
galaxies. 



1 http://www.ucm.es/info/emir/ 



The size of the sources is a relevant parameter in this study, 
given the finite resolution of instruments, and the fact that grav- 
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Fig. 2. Typical error in the magnification factor yu, as a func- 
tion of the magnification for Abell 1835. The solid curve gives 
the statistical error derived from the MCMC model, while the 
dashed curve gives the systematic error between two choices of 
parametrization (see text for details). The thick solid curve is the 
quadratic sum of both errors, used later in the calculation. The 
vertical line represents the conservative upper limit of yu=25. 



itational magnification preserves surface brightness. High red- 
shift galaxies are expec ted to be very small a t z > 7, typically 
< O.'TO on the sky (e.g. Barkana & Loeb 200(3). Recent observa- 
tions of photometric cand idates support this idea dBouwens et alj 
120081 lOeschetalJ 120091) . This means that a large fraction of 
lensed sources should remain spatially unresolved in ground- 
based surveys, even with a strong gravitational magnification 
(hereafter p) of ~ 10. The high resolution capability of JWST 
is clearly needed for resolving such faint galaxies. In the present 
simulations and for detection purposes, we consider all sources 
at z > 6 as spatially unresolved. However, galaxy morphol- 
ogy and image sampling are important when discussing the ef- 
ficiency of surveys based on space facilities, as discussed in 
Sect.lO 



2.1 .2. Lensing effects 

The present simulations address the effect of lensing by a fore- 
ground galaxy cluster. Several well-studied examples of lensing 
clusters are used in order to evaluate the influence of different 
mass distributions on the final results. Reference lensing clusters 
usually display several multiply-imaged systems with redshift 
measurements, allowing us to model their lensing properties ac- 
curately. Lensing clusters considered in these simulations have 
been previously used as gravitational telescopes to search for 
high redshift dropouts and LAEs. We take advantage of this sit- 
uation to perform a direct comparison between our estimates and 
available observations. Finally, we selected clusters with differ- 
ent redshifts, total mass distributions and morphologies, because 
all of these factors are susceptible to affect the way they magnify 
background galaxies. 

We selected three clusters satisfying the previous criteria: 
Abell 1689, Abell 1835 and AO 14. Abell 1689 is one of 
the most spectacular gravitational telescopes, having the largest 
Einstein radius obse rved to date (45"). Bot h optical dropouts 
dBradlev et al.ll2008l) and Lyman-a emitters dStark et akl feOO?) 
candidates have been reported in the background of this clus- 
ter. Abell 1835 and AO 14 are both massive, X-ray luminous 



clusters, previously used in our deep nea r-infrared survey fo r 
high redshift dropouts with VLT/ISAAC dRichard et all 12006). 
Finally, these three clusters constitute the sample used by the 
ZEN2 survey for LAEs in narrow -band images ( Willi s etaTI 
120081) . 

We used the most recent mass models available for the ref- 
erence clusters to derive the magnification maps (see Table Q]) 
although simulation results are found to be weakly sensitive to 
modeling details. Each lensing cluster has been modeled in a 
similar way using the public lensing softwa re Lenstoo fl includ- 
ing the new MCMC optimization method d Julio et all 2007) pro- 
viding bayesian estimates on each parameter derived from the 
model. 

The structure of mass models is given by a sum of indi- 
vidual dark matter subcomponents of two different types: large 
scale components, reproducing the cluster-scale behavior of 
dark matter, and small scale potentials centered on each cluster 
galaxy, reproducing the effect of substructure. Each lensing po- 
tential is parametrized by a Pseudo-Isothermal E lliptical Mass 
Distribution model (PIEMD, [Kassiola & Kovned [19931) . with a 
projected mass density 2 given by: 



r cut 



2G r cut - r core L ( r 2 core + p 2 ) 1 ' 2 (r 2 u , + 



], (2) 



where p 2 = [(jc - jc c )/(1 + e)] 2 + [(y - y c )/(l - e)] 2 , (x c , y c ) 
stands for the central position with respect to the BCG (cD or 
central Bright Cluster Galaxy), e - (a - b)/(a + b) is the ellip- 
ticity, <t () is the central velocity dispersion and (r core ,r cut ) are two 
characteristic radii. For each lensing potential, the position of x 
and y axis is given by the position angle 9. The total mass of 
this profile is proportional to r cu , cr^. The PIEMD parametriza- 
tion, easily linked to the observed geometry of elliptical lensing 
galaxies, has been widely used to model the strong lensing prop- 
erties of massive clu sters (S mith et al.l2 005: Ric hard et al.B2 007: 
Limousin et~alll2007l) . 

A good approximation of the angular distance of the critical 
line, corresponding to maximum magnification in a flat universe, 
is given by the Einstein radius 0e'. 



4?r crl p A { Zs ) - p A ( Zc ) 
c 2 D A (z s ) 



(3) 



where Da(z) stands for the angular distance at redshift z. 
Source and cluster redshifts are, respectively, z s and z c - 

The value of 8e provides a fair estimate of the extension of 
the strongly magnified area (p > 10) in the image plane. This 
value quantifies the power of a gravitational telescope to mag- 
nify background sources. Equation [3] shows that, for a given 
source redshift z s , Qe depends on cr\ and on the cluster redshift 

Zc- 

For the three clusters mentioned before, there is a signifi- 
cant variation in redshift (z c o ~ 0.17 - 0.3) and in <x (taken 
from the mass models and reported in Table[TJ, Abell 1689 being 
~ 30% more massive and less distant than AO 14, for instance. 
We explored a wider range of cluster redshifts in our simulations, 
producing fiducial lensing clusters by adjusting z c between 0.1 
and 0.8 in the three cases, assuming no evolution in the cluster 
properties. This is clearly an over-simplistic and conservative as- 
sumption, as massive clusters of galaxies undergo a dynamical 
evolution during cluster assembly at high redshift. 



http : //www . oamp . fr/cosmology/lenstool 
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Fig. 1. Magnification maps for the three clusters at z s = 8 used in this study (from left to right A1689, A1835 and AC1 14 respec- 
tively). The global size of the image corresponds to the 6' x 6' FOV whereas the blue and red squares correspond respectively 
to 2.2' X 2.2' and 0.85' x 0.85' FOV (see Sect. |2.1.3t . Black contours represent different magnification regimes with increasing 
magnification values from 0.5 to 3 mags towards the cluster center. 



The relevance of the MCMC approach of the cluster mass 
modeling is to derive relevant statistical errors in the magnifica- 
tion factors. Fig. [2] illustrates the typical errors in the magnifi- 
cation at a given position of the lensing field, in the case of the 
lensing cluster Abell 1835. Similar errors in the magnification 
are found in the case of the two other clusters, as all of them 
have > 5 multiple systems constraining independent regions of 
the cluster cores, the majority of them having spectroscopic red- 
shifts. For reasonable magnification factors (< 3 magnitudes), 
this error is always smaller than 0. 1 magnitudes (or ~ 10% rela- 
tive error in flux). For larger magnifications factors, correspond- 
ing to the vicinity (< 2") of the critical lines, the error can reach 
much higher values. The systematic errors in the magnification 
factors, due to the choice of the parametrization when building 
the lensing model, can be estimated for Abell 1835, which have 
been modelled by Richard et al.(2009, submitte d) using both 
PIEM D profiles and Navarro-Frenk- White (NFW. iNavarro et aU 
1 19971) profiles for the cluster-scale mass distributions. The com- 
parison of magnifications from both models, at a given position, 
gives an estimate of the systematic error in the magnification, 
which dominate at large fi (Figure [2]), reaching typical values of 
0.3 magnitudes. We adopted a conservative upper limit of /i = 25 
to avoid singularities in the magnification determination. This is 
justified by the finite resolution of instruments, and the limited 
knowledge on the precise location of the critical lines at such 
high z (typically ~ 1"). The affected area is not significant once 
averaged over the entire field of view. Nevertheless, the quadratic 
sum of the statistical and systematic errors in the magnification 
is later used to derive errors in the number density calculations 
when looking at lensing fields. 

2.1.3. Survey Simulations 

In addition to cluster and source properties, the main ingredients 
to consider in the simulations are the following: 

- The typical field of view (FOV) of near-IR instruments for 
8-10 meters class telescopes and space facilities. The former 
typically range between a few and 10' on a side (e.g. ~ 6' x6' 
for EMIR/GTC, ~ 7'x7' for Hawk-I/VLT). The later are usu- 
ally smaller (e.g. ~ 1' X 1' for NICMOS/HST, ~ 2.2' x 2.2' 
for JWST or WFC3-IR/HST). Fig. Q] presents the compari- 



son between these typical FOV values and the magnification 
regimes found in lensing clusters. The references for the dif- 
ferent instruments used in the simulations are presented in 
Table [TJ 

- The limiting magnitudes of present near-IR surveys based on 
ground-based and space observations tailored to select LBGs 
at z > 6. The former are typically limited to AB ~ 25.5 
(see Sect. 12.1. Il l, whereas the later could reach as deep as 
AB ~ 29 with JWST (see Sect.|33]and |4j. 

The shallow magnitude limit of AB ~ 25.5 achieved on 
ground-based observations should allow us to detect galaxies 
with a UV continuum corresponding to a SFR ~ 4-Q/jj M Q /yr 
at z ~ 10, whereas the typical depth for JWST should be ~ l//i 
M e /yr. 

We can relate the UV luminosities of high redshift galaxies 
with the expected Lyman-ff emission line by convert ing L into 
a star formation rate SFR using the calibrations from iKennicuttl 
(fl998» : 

SFR(M yr _1 ) = 1.4 10~ 28 L (ergs Hz" 1 s" 1 ) (4) 

The expected Lyman-ff luminosity produced at the equilib- 
rium can be written as: 

L tlv (ergs s- 1 ) = (1 - f esc ) f a SFR(M yr" 1 ) (5) 

where f esc is the escape fraction of Lyman-ff photons and 
f a the Lyman-ff production rate per unit of star formation. 
Assuming no reddening, the typic al values for f a range between 
2.44 10 42 and 6.80 10 42 ergs.s~ l dSchaerer 2002). We use these 
scaling relations when discussing the detectability of Lyman-ff 
in lensing fields. 

2.2. Implementation 

We explicitly compute the expected number counts N(z, mo) of 
sources at the redshift z brighter than a limiting magnitude mo 
by a pixel-to-pixel integration of the (magnified) source plane 
as a function of redshift, using the sources and lensing cluster 
models described in the previous subsections. Number counts 
are integrated hereafter within a redshift slice Az = 1 around 
z, unless otherwise indicated. With respect to a blank field, the 
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Parameter 


Explored range 








Reference 


T F 

i^r {a) 


(< z >= 6.0) 


a = 1.6, 


<j>" = 0.4 10~ 3 MpcT 3 , 


M* = -21.07 


DecKwiui et ai. 


(b) 


(< z >= 5.9) 


a = 1.74, 


0* = 1.1 10 3 MpcT 3 , 


AT = -20.24 


riouwens et ai. (2UUs) 


Lb (c) 


(3.8 < z < 7.4) 


a = 1.74, 


<p* = 1.1 10- 3 MpcT 3 , 


M* = -21.02 + 0.36 (z-3.8) 


riouwens et al. (iUUoJ 


Source redshift 


z s = 6 - 12 










Lensing cluster 


Abell 1689 


z c0 =0.178, 


cr =1320 km/s, 


N suh =266 


Limousin et al. (2007) 




Abell 1835 


Z rfD = 0.25, 


cr =1210 km/s, 


Af sub =90 


bmitn(2UUj ) Kicnara(2U0y) 




AC 114 


Zc/o =0.31, 


cr =1080 km/s, 


JV sub =28 


INataraian et ai. {lyya) 


\ 1 LIS Ltil 1 CU-ftlill L 


0.1-0.8 










Survey strategy 


GTC/EMIR 


FOV=6' x 6', 


pix: 0.2" / pixel, 


depth: Az = 1 


Garzon et al. (2007) 




JWST/NIRCam 


FOV=2.2' x 2.2', 


pix: 0.06" 1 pixel, 


depth: Az = 1 


Rieke et al. (2005) 




HST/NICMOS 


FOV=0.85' xO.85'. 


pix: 0.2" / pixel, 


depth: Az = 1 


Thompson (1998) 



Table 1. Summary of the parameters included in o ur simulations. Fo r each entry, we give the range of values explored and reference 
to the relevant publication. LF(a) is the same as in lSteidel et alj dl999l) . but O* is a factor of 3 smaller (see Sect. 12.1. Ti l. 



magnification pushes the limit of integration to fainter magni- 
tudes, whereas the dilution effect reduces the effective volume 
by the same factor. 

An important effect to take into account in cluster fields is 
light contamination coming from the large number of bright 
cluster galaxies, which reduces the surface area reaching the 
maximum depth, and consequently prevents the detection of 
faint objects, especially in the vicinity of the cluster center. This 
contamination effect can be as high as 20% of the total surface 
(Rich ard et al.l l2006h . whereas it is almost negligible in blank 
field surveys. 

We created bright-objects masks by measuring the central 
position (x c ,y c ) and shape parameters (a, b, 6) of galaxies in the 
three cluster fields, each object being appro ximated by an el- 
lipse during this process. We used SExtractor dBertin & Arnoutsl 
1996) in combination with reasonably deep and wide ground- 
based images available from the ESO archive (larger than 6' x 
6' used in these simulations). The characteristics of these im- 
ages are summarized in Table [2] They were reduced using stan- 
dard IRAF routines. The image mask M(x, y) produced is the su- 
perposition of ellipses for all objects in the photometric catalog 
where pixels belonging to object domains were flagged. Ellipses 
correspond to 1 - cr isophotes over the background sky. In other 
words, only images lying on empty regions have been included 
in the lensed samples, thus providing a lower limit for detec- 
tions in lensing fields. This fractional area covered by foreground 
galaxies ranges between 6% and 12% depending on the cluster 
as well as the size and central position of the field of view (Table 
|2]). The largest hidden area corresponds to the smaller field of 
view centered on the cluster (JWST-like). NICMOS pointings 
are even smaller, but they are centered on the critical lines in our 
study and avoid the crowded central regions of the cluster. In 
blank fields, this value doesn't exceed 3-4%. In the next sections, 
we have taken into account this correction in the calculations of 
number counts, both in blank fields and lensing fields. 

Including the object mask M(x,y), number counts N(z,mo) 
are given by the following expression: 



N(z, mo) 



= f J M(x,y) j 



Cv(x, y, z) j L(ju, z, m ) \ a 



(j.(x,y,z) \ L* 

L(p,z,m ) (6) 



exp - 



L(pL, z, mo) \ I L 



where jj(x, y, z) is the magnification induced by the lensing 
field, C v (x,y,z) is the covolume associated with a single spatial 



resolution element pixel with Az = 1 and (L*, (j)*, a ) are the 
parameters of the LF. 

A conservative upper limit of p. — 25 was adopted in the 
vicinity of the critical lines in order to avoid singularities in 
the magnification/dilution determination. This is justified by the 
finite resolution of instruments, and the limited knowledge on 
the precise location of the critical lines at such high-z (typically 
~\"). 

When exploring the impact of cluster redshift, we assumed 
no evolution in the physical parameters r core , r cu t and cr of in- 
dividual potentials, thus keeping the total mass of the cluster 
constant in this process. Variations in z c from the original red- 
shift ZcO of the cluster produce a geometrical effect on the central 
positions (x cj , y cj ) of each PIEMD potential i, measured from a 
reference position (xq, yo) which coincides with the center of the 
BCG: 



XciiZc) -XQ= - — ^ (x ci (Zc ) - X ) 



yaizc) -yo = 



D A (Zc ) 

Da(z c ) 

D A (Zc ) 



(yci(z Co ) - yo) 



(7) 



Similarly, we produced fiducial masks at different cluster 
redshifts z c based on the reference mask at z c o and adjust- 
ing the parameters of the galaxy ellipses, applying the same 
scaling relations on (x c ,y c ). The sizes a and b were scaled by 
D A (z cQ )/D A (z c ). 

3. Results 

As discussed above, lensing introduces two opposite trends on 
the observed sample as compared to blank fields : gravitational 
magnification by a factor p, increasing the expected number of 
sources and thus the total number of objects, and reduction of the 
effective surface by the same factor thus leading to a dilution in 
expec ted counts. This effect was first studied bv lBroadhurst et at] 
(119951) . 

If we consider, for a given redshift z, the cumulative abun- 
dance of sources (per unit of solid angle) with a luminosity 
greater than L and by redshift bin, the magnification bias will 
change depending on p according to 



"L.„„ g (> L,z)=N(> L/p,z)/p(z) 



(8) 
(9) 



where f3 is the logarithmic slope of n(L, z) assuming that this 
function is well represented by a power law in this interval of 
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Cluster Instrument and Reference Filter Exposure Time Depth (5cr, AB) Mask area 6' x 6', 2.2' x 2.2', 0.85' x 0.85' 

Abell 1689 ISAAC/VLT Moorwood (19971 SZ 14.4 ksec 26i0 7% 12% 10% 

Abelll835 FORS/VLT Appenzeller et al. (1998) V 19.6 ksec 27.0 6% 10% 8% 

AO 14 SOFI/NTT Moorwood et al. (1998) K 10.8 ksec 22.0 6% 9% 7% 

Table 2. Characteristics of the images used to produce the foreground object's mask for each cluster field, the three last columns 
representing the fractional surface area covered by the foreground galaxies in each cluster, for the three reference field of view, 
EMIR-like, JWST-like and NICMOS-like respectively 



luminosities : ft = -d(lnn)/d(lnL). The effect on number counts 
is as follows: 

- if /3{z) > 1 the number counts will increase with respect to a 
blank field, and 

- if /3(z) < 1 there will be an opposite trend: i.e. a depletion in 
number counts 

With increasing depth, the /? parameter will decrease in a greater 
or lesser amount depending on the LF, the FOV (because it de- 
termines the mean fi) and the redshift of sources, leading to a 
depletion in number counts in lensing fields as compared to a 
blank fields. With these simple considerations, we expect lens- 
ing clusters to be more efficient than blank fields in relatively 
shallow surveys. 

The efficiency of using lensing clusters as gravitational tele- 
scopes to find high-z galaxies can be quantified with simple as- 
sumptions taking advantage of the properties of the sources ex- 
plained in Sect. |2] In this section, we discuss the results obtained 
by exploring the relevant intervals in the parameter space. We 
present a comparison between the number counts expected in 
lensing and blank fields, as a function of source redshift and for 
different LFs. The influence of lensing cluster properties and red- 
shift is also studied, as well as the expected field to field variance. 

3. 1 . The influence of the field of view 

Here we discuss on the influence of the FOV in the simulations 
for typical surveys. The influence of the limiting magnitude will 
be discussed in Sect. [4] Three different FOV are considered here: 

- 6' x 6' ("EMIR-like" aperture) 

- 2.2' x 2.2' ("JWST-like" or "WFC3/HST" aperture) 

- 0.85' x 0.85' (NICMOS/HST aperture) 

In the last case, the FOV is centered along the critical lines in 
order to achieve the highest mean magnification (Fig.[T). 

The limiting magnitude is AB< 25.5, a value ranging be- 
tween L*(z=6) and 3L« - 5L* at redshift z s ~ 7 to 10. The cluster 
model corresponds to AO 14, but the results are qualitatively 
the same with other models. Fig. [3] displays the relative gain in 
number counts between lensing and blank fields as a function 
of sources redshift, for the three values of the FOV mentioned 
above, and for the three LF adopted in the present simulations. 

The largest gain is obtained for the smallest FOV, as expected 
from geometrical considerations, because the mean magnifica- 
tion decreases with increasing FOV, and in this case fi > 1 given 
the shallow depth (AS; 25.5). For a given FOV, the difference 
between lensing and blank field results strongly depends on the 
shape of the LF. Hence, the comparison between lensing and 
blank field number counts is likely to yield strong constraints on 
the LF, provided that field-to-field variance is small enough. This 
issue is addressed in Sect. 13.51 In the following subsections, we 
adopt a 6' x 6' FOV unless otherwise indicated. 




LF(c) 



°6 8 10 12°6 8 10 12°6 8 10 12 

Z s 

Fig. 3. Relative gain in number counts between lensing and 
blank fields as a function of the source redshift, with Az = 1, 
for different fields of view: 6' x 6' in black, 2.2' x 2.2' in blue 
dotted line and 0.85' x 0.85' in red dot-dashed line (from bottom 
to top). The three panels from left to right represent the 3 values 
of the LF, (a), (b) and (c) respectively. 



3.2. Lensing versus Blank field efficiency 

In this section, we study the effects of lensing clusters on source 
counts, using lensing models for the three reference clusters. 
We compute the expected number of sources brighter than mo, 
the typical apparent magnitude reached in ground-based near-IR 
surveys. The comparison between expected number counts of 
galaxies in a typical 6' x 6' FOV, up to mo < 25.5, per redshift 
bin Az = 1, in a blank field and in the field of a strong lensing 
cluster are presented in Fig.|4]in logarithmic scale. 

We also estimate the error on number counts due to the un- 
certainties on magnification factors (Sect. I2.1.2t . The choice 
of the LF has no influence on the following results. Field to 
field variance dominate the error budget whatever the regime. 
Statistical errors and systematic errors on lensing models are 
smaller but not negligible as their contribution is less sensitive 
than field to field variance to the number of objects. In particu- 
lar, when the number of detected sources is relatively high (i.e. 
when field to field variance is relatively small), they reach ~ 15% 
of the error budget in the worst cases (e.g. for LF(a) at z ~ 6 in 
a 6' x 6' FOV), and typically < 2% when the number of sources 
is relatively small (e.g. for LF(c) at z ~ 8, for any FOV). 

As shown in Fig.|H the presence of a strong lensing cluster 
has a dramatic effect on the observed number of sources, with 
a positive magnification effect. Strong lensing fields are a factor 
between 2 and 10 more efficient than blank fields for the most 
optimistic LF (a), the gain increasing for the LFs (b) and (c), 
reaching a factor between 10 and 100 in the z ~ 6 - 12 domain. 
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A positive magnification bias is observed, increasing with the 
redshift of the sources, and also increasing from optimistic to 
pessimistic values of the LF. This trend is indeed expected given 
the steep shape of the LF around the typical luminosity limits 
achieved in ground-based "shallow" surveys. 

Quantitatively (cf Tableland Fig.|4]l, if the LF for LBGs was 
nearly constant between z ~ 4 and 12, we could always detect 
at least one object over the redshift range of interest. At z ~ 6, 
we expect up to between 7-10 sources, and at z ~ 12 between 
0.7 and 1 galaxies should be detected in a lensing field. Even in 
a blank field, until z ~ 8 - 9 at least one LBG could be found in 
such a large field of view. With more realistic (pessimistic) val- 
ues of the LF (e.g. iBouwens et al.ll2006l 120081) . blank fields are 
particularly inefficient as compared to lensing fields. The size 
of the surveyed area would need to increase by at least a factor 
of ~ 10 in order to reach a number of detections similar to the 
one achieved in a lensing field around z ~ 6 - 8, and this factor 
increases with redshift. Note however that given a limiting (ap- 
parent) magnitude, blank and lensing fields do not explore the 
same intrinsic luminosities (see also Sect|4]i. 

As seen in Fig.|4]and Table[3] there are also some differences 
between the results found in the three lensing clusters, although 
they are smaller than the differences between lensing and blank 
fields for a given LF. The number of expected sources behind 
A 1689 is a factor of two (at z ~ 6) and a factor of three (z ~ 
8) larger than in the other clusters for the realistic LFs (b) and 
(c), whereas the difference is only ~ 10 - 30% for LF (a). The 
influence of lensing properties is studied in Sect. 13.41 

From the results above, it seems that lensing fields allow us 
to detect a larger number of z 1 6 sources based on their UV con- 
tinuum, with some cluster to cluster differences. This result is 
essentially due to the shape of the LF. For magnitude limited 
samples selected within a given field of view, the positive mag- 
nification bias increases with the redshift of the sources and de- 
creases with both the depth of the survey and the size of the sur- 
veyed area. The last trend is purely geometric, as discussed in the 
previous section. The differential behaviour between blank and 
lensing regimes strongly depends on the shape of the LF. The 
comparison between blank and lensing field observations could 
be of potential interest in constraining the LF, provided that field- 
to-field variance is sufficiently small. This issue is addressed in 
the following sections. 

3.3. Redshift of the lensing cluster 

The redshift of the lensing cluster is a potentially important pa- 
rameter when defining an "ideal" sample of gravitational tele- 
scopes. Based on geometrical considerations, we expect the 
magnification bias to decrease with cluster redshift (z c ) after 
reaching a maximum efficiency at some point, depending on 
cluster properties and the size of the surveyed field. The field 
of view considered here is typically a few square arcminutes, 
essentially including the region around the critical lines where 
magnification factors are the highest. Further down, we study 
the impact of z c on the magnification bias. 

Using the non-evolution assumption presented in Sect. 12.21 
we compute the expected number counts for the three refer- 
ence models (A1689, A1835 and AC1 14) with cluster redshifts 
ranging between z=0.1 and 0.8, with a Az = 0.1 step. A step 
Az = 0.05 is used in the z=0. 1-0.3 interval, in order to refine the 
sampling around the maximum value. We use the same depth 
and field size as in previous section. The effect of cluster red- 
shift is clearily seen in Fig. [^representing the number of objects 



as a function of cluster redshift (for the three reference models), 
at a fixed redshift of z=8 for sources. 

The global effect of z c on number counts as a function of the 
source redshift is displayed in Fig. [5] This figure directly com- 
pares to Fig. |4] in the previous section. Table [4] presents the z c 
value which corresponds to a maximum in the expected number 
counts at z=8. This value depends slightly on the source redshift 
and LF. In addition to the Az c = H — 0.05 when changing the 
LF, there is also an increase of z c with higher values of z s , up to 
+0.05 towards z s = 12. The search efficiency of distant galaxies 
in lensing fields is maximised when using clusters at low red- 
shift (z c ~ 0.1 — 0.3). Although the field of view considered here 
is relatively large for near-IR surveys and close to present-day 
cameras, it is the limiting factor at z c ;0.1, where an increasing 
fraction of the strong-magnification area is lost with decreasing 
z c . Also, in this z c regime, the field of view concentrates on the 
central region of the cluster where bright cluster galaxies mask 
an important fraction of the strong-magnification area. The high 
magnification region represents an increasingly small percent- 
age of the field with increasing z c - Number counts in this regime 
asymptotically tend towards a limiting value with increasing z c 
(Fig. IS, wich still is a substantial gain with respect to a blank 
field of the same size. The non-evolution assumption in cluster 
properties has a weak effect on this conclusion. Indeed, clusters 
far from relaxation will be even more ineffective as a gravita- 
tional lenses in the strongest magnification regime. The results 
obtained are hence an optimistic upper limit on number counts 
for realistic clusters beyond the optimal regime z c ~ 0.1 - 0.3. 



l.Or 



LF(a) 



A1689 
A1835 3 
AC114 




LF(b) 

2. =8 



A1689 0.8 
A1835 
AC114 

0.6[ 



0.4 




0.2 



LF(c) 

z, = 8 

A1689 
A1835 
AC114 



0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 0.0 0.2 0.4 0.6 0.8 

Fig. 6. Expected number of objects as a function of the cluster 
redshift for a fixed redshift of sources (z s = 8) and with the 
same depth and field size as in previous section. Three cluster 
models are displayed: A1689 (green dotted line), A1835 (blue 
dot-dashed line) and AC1 14 (red solid line). Panels from left to 
right display respectively LF (a), (b) and (c). 



3.4. Influence of lensing cluster properties 

In this section we focus on the differences between lensing clus- 
ter properties and their influence on expected source counts. As 
seen in previous sections, A1689-like clusters are expected to be 
more efficient irrespective of the cluster redshift. To understand 
this effect, we study the magnification regimes for a reference 
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Table 3. Total number of objects expected within a 6' x 6' FOV (up to Hab < 25.5, Az - 1) from the three LF adopted in these 
simulations. Uncertainties correspond to 1 — cr level in magnification and lensing modeling. 
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Fig. 4. Comparison between the expected number counts of galaxies in a typical 6' x 6' FOV, up to Hab < 25.5, per redshift bin 
Az = 1, in a blank field (dashed lines) and in the field of lensing cluster (solid lines) (left to right respectively A1689, A1835 and 
AC1 14). Expected counts are obtained by the integration of 3 different luminosity functions (a), (b) and (c) from top to bottom. The 
limit of one source detected in the field of view is indicated by an horizontal line to guide the eye. 



Table 4. Redshift of the cluster which maximizes the number of 
objects detected at z=8 for the three LF respectively from top to 
bottom (a), (b) and (c) 



LF 


A1689 


A1835 


AC 114 


actual Zc 


0.184 


0.253 


0.310 


(a) 


0.10 


0.05 


0.25 


(b) 


0.15 


0.10 


0.25 


(c) 


0.20 


0.10 


0.30 



source plane fixed at z s = 8. The distribution of the magnifica- 
tion regimes in the image plane varies from cluster to cluster. 
Histograms in Fig.|7]represent the percentage of the image plane 
(for the 6' x 6' FOV) as a function of the associated magnifi- 
cation. To perform this calculation, cluster redshifts were stan- 
dardized to identical values for a better understanding of the phe- 
nomenon!. As seen in the figure, A 1689 shows a different regime 
at z c ~ 0.8 as compared to the other clusters. While the percent- 
age of the surface affected by strong magnification (p > 10) does 
not exceed 5% in A1835 and AC1 14, it is as high as 8 - 10% in 
A1689, depending on z c - Nevertheless, this difference between 
clusters tends to fade with increasing cluster redshift due to pro- 
jection effects, the fraction of highly magnified pixels becoming 
smaller with respect to the whole FOV. We also note that AC1 14 
and A1835 models have a similar behaviour with minor differ- 



ences (A1835-like clusters being more effective at very low z c 
while the AC1 14 model is more efficient for intermediate z c ). 

Another way of understanding this phenomenom is pre- 
sented in the Fig. [8] where the effective covolume for the 6' X 6' 
FOV is traced as a function of the effective magnitude, for a 
magnitude-limited survey with AB < 25.5. Magnification in 
lensing fields provides an enhanced depth for a magnitude- 
limited survey, where the effective (lensing corrected) covolume 
surveyed decreases with increasing effective depth. The behav- 
ior of A1689 in Fig. [8] illustrates the situation for this particu- 
larly efficient cluster, allowing us to study a ~ lOQMpc 3 vol- 
ume to AB ~ 29.0 with a relatively modest observational in- 
vestment. Except for some particularly efficient lensing clusters 
(such as A 1689), most lensing fields should behave the same 
way as A1835 or AC1 14. 

3.5. Field to field variance 

In this section, we address the expected field-to-field vari- 
ance affecting our previous results in order to estimate its im- 
pact in blank and lensed fields. We used two different ap- 
pr oaches: the two-po i nt corr elation function estimation proposed 
bv lTrenti & Sti avelli (2008) and a pencil beam tracer through the 
Millenium simulation. 

The first estimate i s based on the method implemented by 
iTrenti & Stiavellil (|2008). This method for the calculation of the 
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Fig. 5. The same as Fig. Hbut using different assumptions for the redshift of the clusters, with z c £ [0. 1 , 0.8] and a 0. 1 step (for more 
details see Sect. 13.31 ) 
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Fig. 7. Histogram representing the percentage of the surface (6' x 6' FOV) as a function of the magnification for different redshifts 
of cluster (z c =0.1, 0.2, 0.3, 0.8), using the same color codes for the three clusters as in Fig. |6](A1689: green dotted line, A1835: 
blue dot-dashed line and AC1 14: red solid line) 



cosmic variance is based on the two points correlation function 
£(r) of the sample (Peebles 1993). Field o field variance is given 
by 



where V represents the volume of the survey. 

We focus on the redshift interval z ~ 6 - 8 using the present 
"shallow" survey parameters (see Sect. 14. IK both in blank and 
lensing fields. Here we use the same parameters as in Sect. 13.21 
i.e. typical FOV 6' x 6', m A B = 25.5 and A z = 1. 
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Mabs 
-19.9 -19.4 




Table 6. Number counts for iuab ^ 25.5 and field to field uncer- 
tainties (v f ) calculated from the Millenium simulation in a 6' x6' 
blank field, for different source redshifts 



26.5 27.0 27.5 28.0 

Effective H magnitude. 

Fig. 8. The effective (lensing corrected) covolume sampled at 
z=6.5-7.5 by each cluster is given as a function of effective H AB 
magnitude limit, for a magnitude-limited survey of 6' x 6' FOV 
with Hab < 25.5. The three clusters are displayed with the same 
colors and line codes as in Fig. [6] 



We define the total fract ional error of the counts N follow- 
ing iTrenti & Stiavellil d2008l) (this is the so-called field-to-field 
standard deviation or the, again, improperly called "cosmic vari- 
ance") as : 



v r = 



V<JV 2 > - (AO 2 



(N) 



(11) 



Results are presented in the Table |5]for the three LF consid- 
ered, and for the three typical clusters used in the simulations. 
We note an important field to field variance with such a lim- 
iting magnitude (AB ~ 25.5) either in blank or lensing fields 
due to the small number counts previously derived from calcula- 
tions (see Table |3). Nevertheless, the variance is smaller behind 
gravitational telescopes with the same differential trends men- 
tioned before between the three clusters, i.e. A1689 exhibits a 
stronger magnification bias (see Sect. 13.4-b than the other clus- 
ters which have a similar behavior. Besides, with increasing red- 
shift of sources, the expected number counts decrease leading to 
a larger field to field variance. 

The second estimate was based on the Millennium simula- 
tion, c arried out by the Virg o Co nsortium and described in de- 
tail in Springe fet al.l (120051) and iLemson & Virgo Consortium! 
(2006). The simulation follows N = 2160 3 particles of mass 
8.6 10 s hr x M within a co-moving box of size 500 h 1 Mpc 
on a side. The cosmological model is a ACDM model with 
small differences in the cosmological parameters adopted in 
Sect.Q] but without impact on the final results. These cosmolog- 
ical parameters are consistent with recent determinations from 
the c ombined analysis of the 2dFGRS and 3rd year WMAP 
data (Sanchez & Baugh 2006). Given its high resolution and 
large volume, the Millennium simulation allows us to follow in 
enough details the formation history of a representative sample 
of high redshift galaxy environments. With these prescriptions 
and a realistic beam tracer we can study the field-to-field varia- 
tions in the number counts of star forming galaxies at the epoch 
of interest. 

Our pencil beam tracer is similar to the one developed by 
IKitzbichler & White! d2006l) . We trace through the simulation 
box a parallelepiped where the base is a parallelogram, whose 
size is given by the reference field of view in comoving units, 
and the depth is the comoving depth arbitrarily taken to Az — 1 . 





z s = 6 


z s = 7 


Zs = 8 


(N) 


2.74 


1.10 


0.39 


v r 


68% 


100% 


168% 



The variation of angular distance versus redshift in the redshift 
interval of the selection window considered was properly taken 
into account. This 500/i~' M pc edge box is more than 2000 times 
larger than the effective volume probed by the 6' x 6' FOV : 
~ 10.7 x 10.7 x 323.0(Mpc/h) 3 for instance at z=6. We carried 
out 10000 Monte Carlo realizations of the beam-tracing proce- 
dure by randomly varying the initial position of the beam in or- 
der to calculate the typical number counts of galaxies and the 
associated standard deviation in the field of view with the same 
hypothesis. 

Although this procedure is well suited to determine the field 
to field variance, several studies on this topic suggest an over- 
predi ction on the abundance of m assive galaxies at high redshift 
(e.g. IKitzbichler & White! 120071) . For this reason, we consider 
this second approach as a cross check yielding a lower limit for 
the field to field variance. Results obtained from the Millenium 
simulation are displayed in Table [6] They are in fair agreement 
with those obtained with the first method. 

Field to field variance on number counts obviously depends 
on the depth of the survey. In order to compare our results with 
existing photometric surveys, we calculated the number counts 
of sources in blank and lensing fields (here AC 114) with the 
evolving LF(c) for different deeper magnitude limits (;«/, m = 
27.0, 28.0 and 29.0), in our reference field of view using the 
same parameters as in Sect. 13.21 (see Table |7). The correlation 
function was used to derive the cosmic variance. The total frac- 
tional error (v,) strongly decreases with increasing photometric 
depth, as expected given the increasing number of sources de- 
tected in such a large FOV (e.g. at mu m = 29.0, the total number 
of sources is 1000 times larger than in the "shallow" survey), 
both in blank and lensing fields. The fractional error appears 
slightly larger in lensing than in blank fields at z=6, but this ef- 
fect reverses with increasing source redshift. These estimates for 
the blank field can be compared to present-day surveys. For in- 
stance , the field to field variations obtained by [Bouwens et al.l 
(2006) for a single ACS pointing at z ~ 6 for a limiting magni- 
tude Z85o ~ 29 is 35%. Using the same observational constraints 
(FOV, depth, ...), our simulations yield a v r ~ 30%, a value 
w hich is smaller but fai rly compatible with the results quoted 
bv lBouwens et all (120061) . 

4. Discussion 

4. 1 . Survey parameters and efficiency 



As discussed in Sect. I2.1.3I the FOV and the limiting magni- 
tude are two important survey parameters used in these simula- 
tions. The influence of the FOV for a fixed limiting magnitude 
strongly depends on the shape of the LF. The highest ratio in 
number counts between lensing and blank fields can be achieved 
with the smallest FOV due to simple geometrical considerations. 
This section specifically addresses the evolution on the survey 
efficiency in lensing and blank fields as a function of the limit- 
ing magnitude. 

For these purposes, we use the same approach as in Sect. l3.2l 
to derive number counts within a 6' x 6' FOV in blank fields 
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Table 5. Number counts, field to field variance calculated with the correlation function both in blank and lensing fields, for z = 6, 7 
and 8 within a 6' x6' FOV, for a shallow survey with AB < 25.5. Field to field variance in lensing fields includes 1 - <x magnification 
errors. 
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Table 7. Field to field variance for 3 different magnitude limits : m/,,„ = 27.0, 28.0 and 29.0, in a 6' x 6' blank field and lensing field 
(behind AC1 14) for the LF(c). Field to field variance in lensing fields includes 1 — cr magnification errors. 
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and behind lensing clusters. AC1 14 is used here as a representa- 
tive lensing cluster. Fig.[9]displays the expected number counts 
as a function of the redshift of sources, for different depths 
(AB < 26.0, 27.0, 28.0 and 29.0). An opposite trend between 
blank and lensing fields appears, depending once again on the 
LF and on the redshift of sources. With increasing limiting mag- 
nitude, the efficiency of the survey towards a foreground clus- 
ter diminishes and becomes less effective than in blank fields 
leading to a negative magnification bias for the faintest limiting 
magnitudes (e.g. for LF(a) between AB ~ 21.0 - 28.0, for LF(b) 
between AB ~ 28.0 - 29.0 and for LF(c) beyond AB ~ 29.0). 
This trend, however, is highly sensitive to the FOV. In particular, 
the negative magnification bias appears towards the typical mag- 
nitudes achieved by space facilities (JWST). Fig. [10] displays 
the same results as in Fig. [9] but for a 2.2' x 2.2' FOV (JWST- 
like). The main characteristics remain broadly unchanged, the 
general trends are just exacerbated, the inversion happening to 
lesser depth. 



Lensing and blank field surveys do not explore the same in- 
trinsic luminosities, as shown in Fig. QT| and Q~2] These figures 
compare the expected number density of sources as a function 
of their intrinsic UV luminosity (or equivalent SFR) for differ- 
ent limiting magnitudes ranging from AB;25.5 to 29.0. In the 
case of lensing fields, two different results are given, depending 
on the FOV around the cluster center. In this particular case, the 
source redshift is arbitrarily fixed to z s = 8, assuming a strongly 
evolving LF(c), and the lensing cluster is AC1 14. 

In summary, the number of z>8 sources expected at the typ- 
ical depth of JWST (AB ~ 28 - 29) is much higher in lensing 
than in blank fields if the UV LF is rapidly evolvin g with red- 
shift (LF(c)), as suggested by iBouwens et al.l J2008). The trend 
should be the opposite if the LF remains unchanged between 
z ~ 6 and 8. Lensing clusters are the only way to study the 
faintest building blocks of galaxies, with typical SFR<0.1 to 
\M Q lyr. On the contrary, wide field surveys covering 10 3 to 10 4 
arcmirT 2 are needed to set reliable constraints on the brightest 
part of the LF at z J 6, i.e. for galaxies with SFR > 10M Q /yr. 
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4.2. Influence of galaxy morphology and Image sampling 

Gravitational magnification (e.g. in the tangential direction) in- 
duces an elongation of images along the shear direction while 
preserving the resolution in the perpendicular direction and the 
surface brightness of high redshift galaxies. All the comparisons 
between lensing and blank fields in our simulations assumed that 
observations were conducted with the same intrument setup in 
terms of FOV and spatial sampling, and with the same observa- 
tional conditions, in particular the same limiting surface bright- 
ness and PSF. However, when comparing magnitude-limited 
samples in lensing and blank fields, it is worth discussing the 
influence of galaxy morphology and image sampling on the 
present results. In particular, the evolution in the surface bright- 
ness of high redshift sources is susceptible to hinder the search 
efficiency in clusters if, for instance, number counts in clusters 
were dominated by sources below the limiting surface bright- 
ness. 

As explained in Sect. l2.1.Tl all the previous results have been 
obtained assuming that galaxies at z > 7 are compact as com- 
pared to spatial sampling. Indeed, high redshift sources are ex- 
pected to be very small, typic ally SO.'TO on the sky , based on 
cosmological simulations (e.g. lBarkana & L oeb 2000), in such a 
way that the high resolution capability of JWST is needed for re- 
solving such faint galaxies. Recent observat ions of LBGs candi- 
dates in the HUDF fully support this idea (Bouwens et al. 2008, 
l2009btlOesch et al.ll2009l) . In a recent paper. lOesch et all i200% 
measured the average intrinsic size of z ~ 7-8 LBGs to be 
0.7 + 0.3 kpc. These galaxies are found to be extremely com- 
pact, with very little evolution in their half-light radii between 
z ~ 6 and 7, roughly consistent with galaxies having constant 
comoving sizes, at least within the observed luminosity domain 
~ 0. 1 -lL*(z = 3). Smaller physical sizes are expected for higher 
redshift and/or intrisically fainter galaxies, based on the scal- 
ing of the dark matter halo mass or the disk circular velocity 
dMo et al.lll998l) . This differential trend is actually observed be- 
tween the bright (~ 0.3 - 1L ") and the faint (~ 0.12 - 0.3L*) 
samples of lOesch et al] d2009l) . 

If all high-z galaxies exhibit the same compact and uniform 
morphology, the effective mean surface brightness of a lensed 
galaxy will be brighter or fainter with respect to a blank field 
galaxy with the same apparent magnitude depending on the 
spatial resolution (in practice, the instrumental PSF). The ma- 
jority of lensed sources should remain spatially unresolved on 
their width on seeing-limited ground-based surveys, and even 
on their tangential direction up to a gravitational magnification 
fi ~ 5 - 10. Hence, the apparent surface brightness of a lensed 
source is actually brighter than that of a blank field galaxy of 
similar apparent magnitude (by roughly -2.5log(fi) mags for a 
spatially unresolved galaxy). This situation is typically found in 
the "shallow and wide" near-IR surveys discussed above (e.g. for 
the 6' x6' FOV), where lensing clusters are particularly efficient. 

On the contrary, for a fixed apparent magnitude, the effec- 
tive mean surface brightness of a lensed galaxy is expected to 
become fainter with respect to a blank field galaxy when the 
image resolution is similar or better than its (lensed maximum) 
half-light radius, reaching ~ 2.5log(fi) mags in the worst case. 
This situation is typically expected in the "deep and narrow" 
near-IR surveys with space facilities. In practice, the best spatial 
resolution presently achieved with HST/WFC3 in the near-IR is 
~ O.'TO, reaching 0."065 with JWST/NIRCam, i.e. the typical 
size of the brightest z ~ 7 - 8 LBGs candidates presently iden- 
tified. Therefore, the majority of lensed sources should remain 
spatially unresolved on their width. A lensed source entering the 



apparent-magnitude limited sample because of its magnification 
fi has also a smaller physical size, by a factor of fi (assuming 
a constant M/L scaling with the halo mass) or fi^ 2 (assuming 
a constant M/L scaling with the halo circular velocity), leading 
to an apparent increase on its surface brightness with respect to 
blank-field observations of the same galaxy. Given the spatial 
resolution achieved with HST and JWST, this intrinsic-size ef- 
fect tends to compensate the image dilution described above, in 
such a way that the actual surface brightness of the lensed galaxy 
should get close to the surface brightness of a blank field galaxy 
of similar apparent magnitude. 

For the reasons explained above, and to the best of present 
knowledge, we do not expect the apparent-magnitude limited 
number counts derived in clusters to be strongly biased by 
sources below the limiting surface brightness, provided that the 
usual scalings apply to the size of high-z sources. 

4.3. Comparison with current survey results 

We have compared our simulation results to recent observa- 
tions looking for high-z LBGs. For i nstance, the discover y of a 
bright z = 7.6 + 0.4 lensed galaxy bv lBradlev et al.l d2008l) . with 

AB=24.7 (intrinsic AB 22.6), in a 2.5' x 2.5' FOV survey 

around A1689 is in fair agreement with our expectations. Indeed, 
given the survey characteristics and including ~ 100% variance 
for Zs = 7.6 + 0.4, we expect between 0.2 and 0.8 such bright 
objects in this lensing field, if the LF remains constant between 
z ~ 4 and 8 (LF(a)). In case of a strongly evolving LF(c), the 
expected number of sources in this survey is 0.12 (i.e. ranging 
between and 0.5 with ~ 200% variance) making the discovery 
of this bright source particularly fortunate. Our results for lens- 
ing fields are also cons i stent w ith the number of z ~ 7.5 LBGs 
found by Rich ard et al.l d2008l) . t o the depth of t heir su rvey, us- 
ing LF(b) or (c). Quantitatively, iRichard et al] d2008l) detected 
5 sources with 12 pointings over 6 clusters. With our simula- 
tions, 2.6*qI objects with a variance of ~ 75% are expected with 
the LF(c) model. We also compared with the surface density of 
z ~ 7 ca ndidates in the deep near-IR data behind clust ers ob- 
tained bv lBouwens et all d2009al) . lB~ouwens et all d2009al) found 
a surface density of O.OS^q^ arcmirT 2 with AB<25.5 with a 
typical NICMOS3 FOV. With the strongly evolving LF(c) and 
the same survey characteristics used in our simulations, we ex- 
pect a surface density of 0.01 arcmirT 2 behind a typical cluster 
such as AC114, with a variance of ~ 110%. This result shows 
a relatively good agreement by taking into account field to field 
variance. 



4.4. Lyman Break versus NB searches 

In this section, we discuss on the relative efficiency of blank 
and lensing fields on the detection of LAEs based on either 
NB surveys or the spectroscopic follow up of LBGs at z>6. 
Although the observational effort required to select z > 7 can- 
didates using the dropout technique seems relatively cheap as 
compared to the NB approach, the two approches are comple- 
mentary, as emphasized by the fact that many objects found by 
Lyma n a emission remain weak or undetected in the continuum 
(e.g.|Rhoads& MalhotraU^oH iKodaira et al]|2003l ICubv et all 
l2003llTaniguchi et alj2005l) . A quantitative comparison between 
the properties of LAEs and LBGs at z > 7 within the same 
volume should provide important information on the Lyman a 
transmission, SFR and other properties of these high-z galaxies. 
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Fig. 9. Expected number counts of objects as a function of the 
redshift of sources in a 6' x 6' FOV, for different limiting magni- 
tudes 26.0, 27.0, 28.0 and 29.0 from bottom to top respectively. 
This calculation is provided both in blank (dotted line) and lens- 
ing fields (solid line) (here AC1 14) and for the three LFs (from 
right to left, (a) in red, (b) in blue and (c) in black) 




M AB (1500A) 



Fig. 11. Cumulative surface density of sources as a function of 
their intrinsic UV luminosity, in blank fields (blue solid line) 
and in the lensing fields with FOV ~ 5 arcmin 2 (JWST-like, red 
solid line), for different photometric depths ranging from shal- 
low (AB ~ 25.5) to deep (AB ~ 29.0) surveys and a strongly 
evolving LF(c). The source redshift is arbitrarily fixed at z=8, 
with Az = 1 




0.06 
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Fig. 10. The same as Fig.|9]but for a 2.2' x 2.2' FOV. Some dif- 
ferences appear in comparison with the Fig. [9] For example, the 
total numbers of high-z galaxies expected behind lensing clus- 
ters (solid lines) and the field (dashed lines) are much larger at 
low limiting magnitude (ytiab - 26.0) but this phenomenom is 
reversed for deeper surveys (m AB = 28.0 - 29.0) (see text for 
details). 



Sinc e the pioneering L a rge Area Lyman A l pha S urvey 
(LATA, IRhoads & Malhotral dfoOlh . IRhoads et all (l2003h ). dif- 
ferent NB surveys in blank fields have provided interesting 
galaxy samples in the z~ 5 - 7 interval, e .g. th e large sam- 
ple of Lya emitters at z ~ 5.7 bv lHu et all d2004l). the z= 6.17 



and 6.53 galaxies fo und respectively by ICubv et al.l d2003l) and 
Rhoads et al. (2004), the two z ~ 6.6 galaxies detected by 
Kodaira et al. (2003), and the galaxy at a redshift z=6.96 found 
bv live et alj d2006). In the latter case, which should be repre- 
sentative of z~ 7 samples, the authors used a combination of NB 
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Fig. 12. Cumulative surface density of sources as a function of 
their intrinsic UV luminosity and SFR, in blank fields (blue solid 
line) and in the lensing fields with FOV 6' x 6' (red dotted 
line) and 2.2' x 2.2' (JWST-like, red dashed line), for differ- 
ent photometric depths ranging from shallow (AB ~ 25.5) to 
deep (AB ~ 29.0) surveys and a strongly evolving LF(c). The 
mean magnification over the whole the field is used to derive 
the lensing points, the true distribution is displayed in Fig. QT| 
The source redshift is arbitrarily fixed at z=8, with Az = 1. 
The conversion from absolute magnit ude to SFR is pr ovided in 
Sect. |2.1.3l using the calibrations from Ken nicutildl998l) . 



imaging at 8150A (SuprimeCam) and broad-band photometry in 
the optical bands to select candidates for a subsequent spectro- 
scopic follow up with DEIMOS/Keck. Their confirmation rate 
is relatively high (18 sources out of 26 candidates), leading to 
0.03 sources/arcmin 2 and redshift bin 5z = 0.1. Similar results 
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Fig. 13. Cumulative surface density of observed sources as a 
function of their Lyman alpha luminosity (6' x 6' FOV, A- = 1, 
redshift of sources fixed at 6.6 for the LF(c)). The density is cal- 
culated in blank (blue solid line) and in lensing fields (red solid 
line) for different limiting magnitudes (from right to left: 25.5, 
26.0, 27.0, 28.0 and 29.0). Dashed lines display number counts 
corrected by transmission value e ~ 10% (dashed red line in 
lensing fields and dashed blue line in blank field). For compar- 
ison, raw numbe r counts extracted from the spectroscopic sam- 
ple of LAEs by iKashikawa et al.l (120061) are also given (black 
solid line). The number density derived from live et al.l {2006) 
at z = 6.96 is also indicated, together with corresponding error 
bars (see text for details). As in Fig. Q~2] the magnification used 
to derive the lensing points is averaged over the entire field. 

are reported by Kashika wa et al.l ([2006). All these sources have 
important Lyman a fluxes (a few xl0~ 17 erg crrT 2 s _1 ), and dis- 
play broad Lyman a lines (cr r ~ 200 km/s). A strong evolution is 
found in the numbe r density of LAEs at z > 7 with respect to the 
z~ 5 - 7 interval (live et al.ll2006l IWillis et al.ll2006l ICubv et al.l 
120071) . 

The number of LAEs expected within a sample of LBGs at 
z > 6 can be estimated using the distribution of Lyman-o' equiv- 
alent widt hs derived for the spectroscopic sample of LBGs at 
Z ~ 3 by Sha plev et al.l (120031 assuming no evolution in the 
population of LAEs with respect to LBGs. This simplistic scal- 
ing should be enough for the simulation needs. We introduce 
a factor e, defined below, which can be linked to the Lyman a 
transmission as follows: 



L Lya (l + Z) W L ya 

where L1500 is the UV monochromatic luminosity at 1500 A, 
LLya is the Lyman-ff luminosity and W]_ ya is the Lyman-ff equiv- 
alent width. With this simple assumption, the average value 
for the Lyman-ff equivalent width is ~ 10 A, corresponding to 
e ~ 10%. This value can be used to derive a rough estimate of 
expected number density of LAEs, from a population of LBGs 
In addition, the number density is also corrected to take into 
account of the fraction of the LBGs sample displaying Lyff in 
emission. 



Fig. Qj] displays the cumulative number counts of sources 
at z ~ 6.6 integrated from the LF(c) as a function of the 
Lyman-ff luminosity, scaled according to the UV luminosities 
(cf Sect. 12. Lit in the typical 6' x 6' FOV, together with a com- 
parison of observations in a similar redshift dom ain (z s ~ 6.6 f or 
Kashikawa spectroscopic sample of LAEs and live et a l. 2006 at 
z=6.96). 

The number density of LBGs at z=7 (with 6z = 0.1, 
close to the band-width of NB surveys) ranges between 0.001 
(LF(c)) and 0.02 (LF(a)) sources/orcm/n 2 for a survey limited 
to H(AB) $25.5, depending on the LF. Lensing clusters improve 
these numbers by a factor ranging between 6 (for LF(c)) and 2 
(for LF(a)). In case of a deep survey limited to H(AB) < 29.0, the 
number densities reach 1 (LF(c)) to 2 (LF(a)) sources/arcm/n 2 . 
In this case, there is a negative magnification bias of the or- 
der of 20%. These numbers, obtained with a simplistic model, 
are between a factor of ~ 10 (for bright sources) and a few 
(fo r faint sources) smaller than the number densities obtained 
by IKashikawa et al.l (120061) for their spectroscopic sample. With 
increasing z s (see Fig. |9j for instance at z=9 with the strongly 
evolving LF(c), no sources can be detected for a shallow sur- 
vey limited to H(AB)< 25. 5 and for a deeper limited survey 
(H{AB) < 29.0), a minimum of 3 arcmin 2 surveyed area is needed 
to obtain 1 source in a blank field. In a lensing field with the 
(LF(c)), these number densities reach 0.002 for H(AB) <25.5 
and 0.32 sources/arcm/n 2 for H(AB)< 29.0. The relatively low- 
efficiency of lensing clusters with NB techniques in t he z > 9 do- 
main has been recently confirmed by the results of Will is et all 
d2008h . 

The preselection of z > 6 - 7 candidates in lensing fields has 
two main advantages with respect to blank fields. In the shallow 
(AB < 25.5) regime, there is an increase by a factor ~ 8 - 10 on 
the number of sources detected and a moderate gain in depth for 
a given exposure time (i.e. ~ 0.5 magnitudes at AB ~ 25.5). In 
the deep-survey regime (AB < 28-29), there is a gain in intrinsic 
depth, for a number of candidates which remains essentially con- 
stant (i.e. ~ 0.8 mags gain at AB < 28). The relative efficiency 
of lensing with respect to blank field counts in Fig. [T2l depends 
on the FOV. The two predictions get close to each other with 
increasing values of the FOV in lensing surveys, and the trend 
goes in the opposite direction for smaller FOV. This trend is the 
same for both LBGs and LAEs. To explore the bright end of the 
LF, blank field surveys are needed with a large FOV, whereas 
lensing clusters are particularly useful to explore the faint end of 
the LF. This trend is further discussed in the next Section. 

4.5. Towards the ideal survey: constraining the Luminosity 
Function of high-z sources 

All present photometric surveys aimed at constraining the UV 
LF at z>7, either space or ground-based, are still dramatically 
small in terms of effective surface. Wide and deep optical+ 
near-IR surveys in lensing and blank fields are needed to set 
strong constraints on the LF and on the star-formation density 
at z>7. An important issue is the combination between photo- 
metric depth and surveyed area which is needed to identify a 
representative number of photometric candidates, or to reach a 
significant non-detection limit in order to constrain the LF of 
z>7 sources. 

There are three different aspects to consider when design- 
ing an "ideal" survey aiming at constraining the LF: the depth 
and the area of the survey, and the corresponding field to field 
variance. In order to address these issues, we have computed 
the expected field to field variance corresponding to lensing and 
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blank field surveys, for different survey configurations (area and 
depth). A summary of these results is given in Table[8]for differ- 
ent number of lensing clusters, and for two representative depths 
in the H-band (i.e. a "shallow" survey with AB < 25.5, and a 
"deep" survey with AB < 28.0) assuming a strongly evolving 
LF(c) in all cases. This table complements the results given in 
Tables[5]and[7]for blank and lensing fields as a function of depth. 
In all cases, we use AC1 14 as a reference for lensing clusters. 

Regarding field-to-field variance in number counts, results 
are expected to be similar in blank and lensing fields for a rela- 
tively wide FOV (~ 40 - 50 arcmin 2 ; see Sect. l3.5l and Table|7]l. 
As shown in Table [8] a deep lensing survey using ~ 10 clusters 
should be able to reach a variance < 20% on sources at 6 < z < 8, 
irrespective of the actual LF. This value is better than present-day 
photometric survey s in blank fields, typ ically reaching 30-35% 
for AB $29.0 (e.g. lBouwens et aLll2008l) . which in turn is rather 
close to what could be achieved in a single lensing cluster for 
AB <28.0. 

A different survey strategy consists of increasing the number 
of lensing clusters with a shallow limiting magnitude. In this 
case, a few tens of lensing clusters (typically between 10 and 
50, depending on the LF) are needed to reach a variance of ~ 
30% at z 5 8. Note that the difference in exposure time between 
the shallow and deep surveys reported in Table [8] is a factor ~ 
100, and that ~ 10 pointings are needed on blank fields in order 
to reach the same number of z ~ 6 - 8 sources as in a single 
"shallow" lensing field. 

In the case of a strongly evolving LF(c), photometric surveys 
should reach a minimum depth of AB ~ 28 to achieve fair statis- 
tics on z ~ 9-10 sources using a lensing cluster (Table|7]i. In this 
case we expect between 20 (z=9, ~ 40% variance) and 8 (z=10, 
~ 40% variance) sources per lensing cluster in a ~ 40 arcmin 2 
FOV. The efficiency is a factor of 10 smaller at z ~ 10 in blank 
fields. Fair statistics at z ~ 12 should require a minimum depth 
close to AB ~ 29 both in lensing and blank fields. 

Constraining the LF of star forming galaxies at z > 7 should 
require the combination of blank and lensing field observations. 
This is illustrated for instance in Fig.fTTIand [T2]for an example 
at z s = 8. A survey towards a lensing cluster has several advan- 
tages. It increases the total number of sources available for spec- 
troscopic follow up, and it helps extending the sample towards 
the faint edge of the LF and towards the highest possible limits 
in redshift. On the other hand, blank fields are needed to achieve 
fair statistics on the bright edge of the LF. Thus an "ideal" survey 
should combine both blank and lensing fields. Given the num- 
bers presented in previous sections, a blank field used for these 
purposes should be a factor ranging between -10 and 100 times 
larger than a lensing field (depending on the redshift domain, 
photometric depth, and actual LF) in order to efficiently com- 
plete the survey towards L > L*. This should be possible with 
the new upcoming surveys, such as the WIRCAM ultra deep 
survey (WUDS) at CFHT (~ 400 arcmin 2 , with YJHK <25.5), 
UKIDSS-UDS (~ 2700 arcmin 2 , with YJHK< 25) or Ultra- 
Vista (~ 2600 arcmin 2 , with YJH< 26, K<25.6). The optimum 
number of lensing fields ranges between ~ 10-20 (for z ~ 6 - 8 
studies with "shallow" photometry) to a few (for "deep" surveys 
targeting z ~ 8 - 12 sources). 



5. Summary and conclusions 

We have evaluated the relative efficiency of lensing clusters with 
respect to blank fields in the identification and study of z > 6 
galaxies. The main conclusions of this study are given below. 



For magnitude-limited samples of LBGs at z > 6, the magni- 
fication bias increases with the redshift of sources and decreases 
with both the depth of the survey and the size of the surveyed 
area. Given the typical near-IR FOV in lensing fields, the max- 
imum efficiency is reached for clusters at z ~ 0.1 - 0.3, with 
maximum cluster-to-cluster differences ranging between 30 and 
50% in number counts, depending on the redshift of sources and 
the LF. 

The relative efficiency of lensing with respect to blank fields 
strongly depends on the shape of the LF, for a given photometric 
depth and FOV. The comparison between lensing and blank field 
number counts is likely to yield strong constraints on the LF. 

The presence of a strong-lensing cluster along the line of 
sight has a dramatic effect on the observed number of sources, 
with a positive magnification effect in typical ground-based 
"shallow" surveys (AB < 25.5). The postive magnification bias 
increases with the redshift of sources, and also from optimistic to 
pessimistic values of the LF. In case of a str ongly evolving LF at 
z > 7, as proposed by Bo uwens et al.l d2008). blank fields are par- 
ticularly inefficient as compared to lensing fields. For instance, 
the size of the surveyed area in ground-based observations would 
need to increase by a factor of ~ 10 in blank fields with respect 
to a typical ~ 30-40 arcmin 2 survey in a lensing field, in or- 
der to reach the same number of detections at z ~ 6-8, and 
this merit factor increases with redshift. All these results have 
been obtained assuming that number counts derived in clusters 
are not dominated by sources below the limiting surface bright- 
ness of observations, which in turn depends on the reliability of 
the usual scalings applied to the size of high-z sources. 

Ground-based "shallow" surveys are dominated by field-to- 
field variance reaching ~ 30 to 50% in number counts between 
z~6 and 8 in a unique ~ 30-40 arcmin 2 lensing field survey (or 
in a 400 arcmin 2 blank field), assuming a strongly evolving LF. 

The number of z>8 sources expected at the typical depth of 
JWST (AB ~ 28 - 29) is much higher in lensing than in blank 
fields if the UV LF is rapidly evolving with redshift (i.e. a factor 
of ~ 10 at z ~ 10 with AB < 28). 

Blank field surveys with a large FOV are needed to probe the 
bright edge of the LF at z > 6-7, whereas lensing clusters are 
particularly useful to explore the mid to faint end of the LF. 
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ABSTRACT 

Context. The identification and study of the first galaxies remains one of the most exciting topics in observational cosmology. The 
determination of the best possible observing strategies is a very important choice in order to build up a representative sample of 
spectroscopically confirmed sources at high-z (z > 7), beyond the limits of present-day observations. 

Aims. This paper is intended to precisely adress the relative efficiency of lensing and blank fields in the identification and study of 
galaxies at 6; z; 12. 

Methods. The detection efficiency and field-to-field variance are estimated from direct simulations of both blank and lensing fields 
observations. Present known luminosity functions in the UV are used to determine the expected distribution and properties of distant 
samples at z > 6 for a variety of survey configurations. Different models for well known lensing clusters are used to simulate in details 
the magnification and dilution effects on the backgound distant population of galaxies. 

Results. The presence of a strong-lensing cluster along the line of sight has a dramatic effect on the number of observed sources, 
with a positive magnification bias in typical ground-based "shallow" surveys (AS; 25.5). The positive magnification bias increases 
with the redshift of sources and decreases with both depth of the survey and the size of the surveyed area. The maximum efficiency 
is reached for lensing clusters at z ~ 0.1 - 0.3. Observing blank fields in shallow surveys is particularly inefficient as compared to 
lensing fields if the UV LF for LBGs is strongly evolving at z i 7. Also in this case, the number of z > 8 sources expected at the typical 
depth of JWST (AB ~ 28 - 29) is much higher in lensing than in blank fields (e.g. a factor of ~ 10 for AS; 28). All these results 
have been obtained assuming that number counts derived in clusters are not dominated by sources below the limiting surface 
brightness of observations, which in turn depends on the reliability of the usual scalings applied to the size of high-z sources. 
Conclusions. Blank field surveys with a large field of view are needed to prove the bright end of the LF at z i 6 — 7, whereas lensing 
clusters are particularly useful for exploring the mid to faint end of the LF. 

Key words, galaxies : formation - galaxies : high redshift - galaxies : photometry - galaxies : clusters : lensing - 



i 1. Introduction high-z up to the limits of their survey, in particular towards the 

J\ bright end of the Luminosity Function (LF). A strong evolution 

• Constraining the abundance of z > 7 sources remains an impor- is alsQ observed m the number densit of sources detected with 

Cd ; tant challenge of modern cosmology. Recent WMAP results place NB techni which seems to be much smaller aU > 7 than in 

the first building blocks of galaxies at redshifts z = 11.0 ± 1.4, ^ z ^ ^ _ 7 m t erva i (???) 

suggesting that reionization was an extended process (?). 

Distant star-forming systems could have been responsible for Both dropout and NB approaches require a subsequent spec- 
a significant part of the cosmic reionization. Considerable ad- troscopic confirmation of the selected candidates. For now ap- 
vances have been made during the last years in the observation of proximately ten galaxies beyond z~6.5 are known with se- 
the early Universe with the discovery of galaxies at -6-7, close cure spectroscopic redshifts (?????). All samples beyond this 
to the end of the reionization epoch (e.g. ??????????), and the redshift are mainly supported by photometric considerations 
first prospects up to z~10 (??????). (??????). This situation is expected to dramatically improve in 
High-z surveys are mainly based on two different and com- the near future with the arrival of a new generation of multi- 
plementary techniques: the dropout (Lyman-a* Break) identifica- ob J ect spectrographs in the near-IR, such as MOIRCS/Subaru, 
tion, which is an extrapolation of the drop-out technique used for Flamingos2/Gemim-S (~ 2009), or EMIR/GTQJ (~ 2012), with 
Lyman-Break Galaxies (LBGs, ?) to higher redshifts (e.g. ???), wel1 sulted field of Vlew ' s P ectral resolution and sensitivity, 
and the narrow-band (NB) imaging aimed at detecting Lyman These forthcoming facilities should provide spectroscopic Con- 
or emitters (LAEs, e.g. ?????). Using the former technique, ? nrmation for a large number of z>7 candidates identified from 
found a strong evolution in the abundance of galaxies between dee P photometric surveys, and the first characterization of the 
z~7-8 and z~3-4, the SFR density beeing much smaller at very 
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physical properties of these sources (e.g. IMF, stellar popula- 
tions, fraction of AGN, ...). 

The aim of this paper is to determine the best possible ob- 
serving strategies in order to build up a representative sample 
of spectroscopically confirmed z > 7 galaxies. The photomet- 
ric pre-selection of high-z candidates could be achieved either 
in blank fields or in lensing clusters. This later technique, also 
first referred to as the "gravitational telescope" by Zwicky , has 
proven highly successful in identifying a large fraction of the 
most distant galaxies known today thanks to magnifications by 
typically 1-3 magnitudes in the cluster core (e.g. ??????). The 
presence of a strong lensing cluster in the surveyed field intro- 
duces two opposite effects on number counts as compared to 
blank fields. In one hand, gravitational magnification increases 
the number of faint sources by improving the detection towards 
the faint end of the LF. On the other hand, the reduction of the 
effective surface by the same magnification factor leads to a dilu- 
tion in observed counts. The global positive/negative magnifica- 
tion bias obviously depends on the slope of the number counts, 
as pointed out by ?. 

This paper addresses the relative efficiency of surveys con- 
ducted on blank and lensing fields as a function of the relevant 
parameters, namely the redshift of the sources, the redshift and 
properties of the lensing clusters and the survey characteristics 
(i.e. area, photometric depth...). This calculation requires a de- 
tailed simulation of observations using lensing models, and re- 
alistic assumptions for the properties of background sources ac- 
cording to present-day observational results, in particular for the 
luminosity function and typical sizes of z > 7 galaxies. 

The paper is organized as follows. In Section|2]we describe 
the simulations performed in order to determine the relative de- 
tection efficiency for high-z sources, both in lensing and blank 
fields. Section [3] presents the results, in particular the detection 
efficiency achieved as a function of redshift for both sources and 
lensing clusters, together with a discussion on the influence of 
lensing cluster properties and field-to-field variance. A discus- 
sion is presented in Section|4]on the relative efficiency as a func- 
tion of survey parameters, and a comparison between simula- 
tions and present surveys. Conclusions are given in Section[5] 

Throughout this paper, we adopt a concordance cosmo- 
logical model, with £2a = 0.7, Q m = 0.3, and Hq 
70 km s Mpc~ l . All magnitudes are given in the AB system. 
Conversion values between Vega and AB systems for the filters 
used in this paper are typically Cab- 0.95, 1.41 and 1.87 in /, H 
and K respectively, with m AB = m Ve ga + Cab- 



2. Simulations of lensing and blank field 
observations 

2.1. Simulation parameters 

This section describes the ingredients used in the simulations to 
implement different assumptions that would affect our efficiency 
in detecting high redshift galaxies. There are three important as- 
pects to be considered in the comparison between lensing and 
blank fields. The first one is the LF and typical sizes of sources. 
The second one concerns the properties of the lensing clusters, 
in particular their mass distribution and redshift. The third one is 
related to the survey parameters, namely the photometric depth 
and the size of the field. All these aspects are discussed in this 
section. Table[T]provides the list of parameters used in these sim- 
ulations, together with the range of values explored. 



2.1.1. Source Properties 

These simulations are focused on the detection of sources in the 
redshift range 6 < z < 12, a relevant domain for spectroscopic 
follow-up with near-infrared instruments. The lower limit of this 
redshift domain overlaps with current photometric surveys mea- 
suring the LF at z ~ 6-7 (e.g. ?). However, the LF is still 
largely unconstrained beyond z S 7 because of the lack of spec- 
troscopic confirmation of photometric samples and the relatively 
small size of the surveyed volumes. 

The abundance of background sources at these redshifts is 
given by the luminosity function 4>(L), with L the rest-frame UV 
luminosity at 1500 A. (p(L) is the most basic description of the 
galaxy population from an observer point of view. We adopt a 
parametrization based on the analytical Schechter function (?): 

<f,(L)dL = f(^Y ! exp (-A)^) (i) 

The slope at faint luminosities a, the characteristic luminos- 
ity L* and the normalization factor O* have been constrained 
by several photometric surveys targeting Lyman Break Galaxies 
(LBG) at high redshift (z > 4) (?, ?, ?, ?, ?, ?). Three different 
representative cases are being discussed in our simulations, ba- 
sically exploring our present knowledge (or lack of knowledge) 
of the overall shape of the LF at z > 6 (Table [TJ: 

(a) An "optimistic" scenario where LBGs show no-evolution 
from z ~ 6, with the LF parameters as determined by ?. 
Indeed, the LF at z ~ 6 found by these authors display the 
same shape as for z ~ 3 (?), but a 3 times smaller normaliza- 
tion factor (but see, e.g., ?). 

(b) A constant LF based on robust measurements by ? at z ~ 6 
in the Hubble UDF, but using the more recent fit parameters 
from ?. As compared to model (a), this LF exhibits a turnover 
towards the bright end. 

(c) The evolutionary LF recently proposed by ?, which includes 
an important dimming of L* with increasing redshift. This 
LF represents the "pessimistic" case with respect to the case 
(a), with very few high-luminosity galaxies. 

The size of the sources is a relevant parameter in this study, 
given the finite resolution of instruments, and the fact that grav- 
itational magnification preserves surface brightness. High red- 
shift galaxies are expected to be very small at z > 7, typically 
<0.'T0 on the sky (e.g. ?). Recent observations of photomet- 
ric candidates support this idea (?, ?). This means that a large 
fraction of lensed sources should remain spatially unresolved in 
ground-based surveys, even with a strong gravitational magnifi- 
cation (hereafter jj.) of ~ 10. The high resolution capability of 
JWST is clearly needed for resolving such faint galaxies. In the 
present simulations and for detection purposes, we consider all 
sources at z > 6 as spatially unresolved. However, galaxy mor- 
phology and image sampling are important when discussing 
the effi cienc y of surveys based on space facilities, as discussed 
in Sect. 1431 

2.1 .2. Lensing effects 

The present simulations address the effect of lensing by a fore- 
ground galaxy cluster. Several well-studied examples of lensing 
clusters are used in order to evaluate the influence of different 
mass distributions on the final results. Reference lensing clusters 
usually display several multiply-imaged systems with redshift 
measurements, allowing us to model their lensing properties ac- 
curately. Lensing clusters considered in these simulations have 
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Fig. 1. Magnification maps for the three clusters at z s = 8 used in this study (from left to right A1689, A1835 and AC1 14 respec- 
tively). The global size of the image corresponds to the 6' x 6' FOV whereas the blue and red squares correspond respectively 
to 2.2' x 2.2' and 0.85' x 0.85' FOV (see Sect. 12.1.3V Black contours represent different magnification regimes with increasing 
magnification values from 0.5 to 3 mags towards the cluster center. 
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Fig. 2. Typical error in the magnification factor /i, as a func- 
tion of the magnification for Abell 1835. The solid curve gives 
the statistical error derived from the MCMC model, while the 
dashed curve gives the systematic error between two choices of 
parametrization (see text for details). The thick solid curve is the 
quadratic sum of both errors, used later in the calculation. The 
vertical line represents the conservative upper limit of yu=25. 



been previously used as gravitational telescopes to search for 
high redshift dropouts and LAEs. We take advantage of this sit- 
uation to perform a direct comparison between our estimates and 
available observations. Finally, we selected clusters with differ- 
ent redshifts, total mass distributions and morphologies, because 
all of these factors are susceptible to affect the way they magnify 
background galaxies. 

We selected three clusters satisfying the previous criteria: 
Abell 1689, Abell 1835 and AO 14. Abell 1689 is one of 
the most spectacular gravitational telescopes, having the largest 
Einstein radius observed to date (45"). Both optical dropouts (?) 
and Lyman-a emitters (?) candidates have been reported in the 
background of this cluster. Abell 1835 and AC1 14 are both mas- 
sive, X-ray luminous clusters, previously used in our deep near- 
infrared survey for high redshift dropouts with VET/ISAAC (?). 



Finally, these three clusters constitute the sample used by the 
ZEN2 survey for LAEs in narrow-band images (?). 

We used the most recent mass models available for the ref- 
erence clusters to derive the magnification maps (see Table Q]) 
although simulation results are found to be weakly sensitive to 
modeling details. Each lensing cluster has been modeled in a 
similar way using the public lensing software LenstooQ includ- 
ing the new MCMC optimization method (?) providing bayesian 
estimates on each parameter derived from the model. 

The structure of mass models is given by a sum of indi- 
vidual dark matter subcomponents of two different types: large 
scale components, reproducing the cluster-scale behavior of 
dark matter, and small scale potentials centered on each clus- 
ter galaxy, reproducing the effect of substructure. Each lensing 
potential is parametrized by a Pseudo-Isothermal Elliptical Mass 
Distribution model (PIEMD, ?), with a projected mass density 2 
given by: 



Z(x,y) = 



r cut 



2G 



(d„+p 2 ) 1/2 (d, + P 2 ) 1/2 



ml (2) 



where p 2 = [(jc - jc c )/(1 + e)] 2 + [(y - y c )l{\ - e)] 2 , (x c , y c ) 
stands for the central position with respect to the BCG (cD or 
central Bright Cluster Galaxy), e - (a - b)/(a + b) is the ellip- 
ticity, <To is the central velocity dispersion and {r core ,r cut ) are two 
characteristic radii. For each lensing potential, the position of x 
and y axis is given by the position angle 8. The total mass of 
this profile is proportional to r cut cr 2 . The PIEMD parametriza- 
tion, easily linked to the observed geometry of elliptical lensing 
galaxies, has been widely used to model the strong lensing prop- 
erties of massive clusters (???). 

A good approximation of the angular distance of the critical 
line, corresponding to maximum magnification in a flat universe, 
is given by the Einstein radius 0e'- 



4tt o% D A { Zs ) - D A (zc) 



(3) 



c 2 D A (z s ) 

where D A {z) stands for the angular distance at redshift z- 
Source and cluster redshifts are, respectively, z s and z c - 



http : //www . oamp . fr/cosmology/lenstool 
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The value of Qe provides a fair estimate of the extension of 
the strongly magnified area (// > 10) in the image plane. This 
value quantifies the power of a gravitational telescope to mag- 
nify background sources. Equation [3] shows that, for a given 
source redshift z s , 8e depends on cr^ and on the cluster redshift 

Zc- 

For the three clusters mentioned before, there is a signifi- 
cant variation in redshift (z c o ~ 0.17 - 0.3) and in cr (taken 
from the mass models and reported in TableQ]), Abell 1689 being 
~ 30% more massive and less distant than AC1 14, for instance. 
We explored a wider range of cluster redshifts in our simulations, 
producing fiducial lensing clusters by adjusting z c between 0. 1 
and 0.8 in the three cases, assuming no evolution in the cluster 
properties. This is clearly an over-simplistic and conservative as- 
sumption, as massive clusters of galaxies undergo a dynamical 
evolution during cluster assembly at high redshift. 

The relevance of the MCMC approach of the cluster mass 
modeling is to derive relevant statistical errors in the magnifica- 
tion factors. Fig. [2] illustrates the typical errors in the magnifi- 
cation at a given position of the lensing field, in the case of the 
lensing cluster Abell 1835. Similar errors in the magnification 
are found in the case of the two other clusters, as all of them 
have > 5 multiple systems constraining independent regions of 
the cluster cores, the majority of them having spectroscopic red- 
shifts. For reasonable magnification factors (< 3 magnitudes), 
this error is always smaller than 0.1 magnitudes (or ~ 10% rela- 
tive error in flux). For larger magnifications factors, correspond- 
ing to the vicinity (< 2") of the critical lines, the error can reach 
much higher values. The systematic errors in the magnification 
factors, due to the choice of the parametrization when building 
the lensing model, can be estimated for Abell 1835, which have 
been modelled by Richard et al.(2009, submitted) using both 
PIEMD profiles and Navarro-Frenk- White (NFW, ?) profiles for 
the cluster-scale mass distributions. The comparison of magnifi- 
cations from both models, at a given position, gives an estimate 
of the systematic error in the magnification, which dominate at 
large y. (Figure|2]), reaching typical values of 0.3 magnitudes. We 
adopted a conservative upper limit of /u — 25 to avoid singular- 
ities in the magnification determination. This is justified by the 
finite resolution of instruments, and the limited knowledge on 
the precise location of the critical lines at such high z (typically 
~ 1")- The affected area is not significant once averaged over 
the entire field of view. Nevertheless, the quadratic sum of the 
statistical and systematic errors in the magnification is later used 
to derive errors in the number density calculations when looking 
at lensing fields. 

2.1.3. Survey Simulations 

In addition to cluster and source properties, the main ingredients 
to consider in the simulations are the following: 

- The typical field of view (FOV) of near-IR instruments for 
8-10 meters class telescopes and space facilities. The former 
typically range between a few and 10' on a side (e.g. ~ 6' x6' 
for EMIR/GTC, ~ 7'x7' for Hawk-I/VLT). The later are usu- 
ally smaller (e.g. ~ 1' x 1' for NICMOS/HST, ~ 2.2' x 2.2' 
for JWST or WFC3-IR/HST). Fig. Q] presents the compari- 
son between these typical FOV values and the magnification 
regimes found in lensing clusters. The references for the dif- 
ferent instruments used in the simulations are presented in 
Table [Tj 

- The limiting magnitudes of present near-IR surveys based on 
ground-based and space observations tailored to select LBGs 



at z > 6. The former are typically limited to AB ~ 25.5 
(see Sect. 12.1. IK whereas the later could reach as deep as 
AB ~ 29 with JWST (see Sect. [33]and 0. 

The shallow magnitude limit of AB ~ 25.5 achieved on 
ground-based observations should allow us to detect galaxies 
with a UV continuum corresponding to a SFR ~ 40/// M Q /yr 
at z ~ 10, whereas the typical depth for JWST should be ~ 1/yiz 
M Q /yr. 

We can relate the UV luminosities of high redshift galaxies 
with the expected Lyman-ff emission line by converting L into a 
star formation rate SFR using the calibrations from ?: 

SFR(M yr 1 ) = 1.4 10~ 28 L (ergs Hz" 1 s" 1 ) (4) 

The expected Lyman-ff luminosity produced at the equilib- 
rium can be written as: 

L iv „(ergs s- 1 ) = (1 - f esc ) f a SFR(M yr 1 ) (5) 

where f esc is the escape fraction of Lyman-o- photons and 
f a the Lyman-a production rate per unit of star formation. 
Assuming no reddening, the typical values for f a range between 
2.44 10 42 and 6.80 10 42 ergs.s' 1 (?). We use these scaling re- 
lations when discussing the detectability of Lyman-Q' in lensing 
fields. 

2.2. Implementation 

We explicitly compute the expected number counts N(z, trio) of 
sources at the redshift z brighter than a limiting magnitude mo by 
a pixel-to-pixel integration of the (magnified) source plane as a 
function of redshift, using the sources and lensing cluster mod- 
els described in the previous subsections. Number counts are 
integrated hereafter within a redshift slice Az — 1 around z, 
unless otherwise indicated. With respect to a blank field, the 
magnification pushes the limit of integration to fainter magni- 
tudes, whereas the dilution effect reduces the effective volume 
by the same factor. 

An important effect to take into account in cluster fields is 
light contamination coming from the large number of bright 
cluster galaxies, which reduces the surface area reaching the 
maximum depth, and consequently prevents the detection of 
faint objects, especially in the vicinity of the cluster center. This 
contamination effect can be as high as 20% of the total surface 
(?), whereas it is almost negligible in blank field surveys. 

We created bright-objects masks by measuring the central 
position (x c ,y c ) and shape parameters (a, b, 0) of galaxies in the 
three cluster fields, each object being approximated by an ellipse 
during this process. We used SExtractor (?) in combination with 
reasonably deep and wide ground-based images available from 
the ESO archive (larger than 6' x 6' used in these simulations). 
The characteristics of these images are summarized in Table [2] 
They were reduced using standard IRAF routines. The image 
mask M(x, y) produced is the superposition of ellipses for all ob- 
jects in the photometric catalog where pixels belonging to object 
domains were flagged. Ellipses correspond to 1 - <x isophotes 
over the background sky. In other words, only images lying on 
empty regions have been included in the lensed samples, thus 
providing a lower limit for detections in lensing fields. This frac- 
tional area covered by foreground galaxies ranges between 6% 
and 12% depending on the cluster as well as the size and cen- 
tral position of the field of view (Table |2|. The largest hidden 
area corresponds to the smaller field of view centered on the 
cluster (JWST-like). NICMOS pointings are even smaller, but 
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Table 1. Summary of the parameters included in our simulations. For each entry, we give the range of values explored and reference 
to the relevant publication. LF(a) is the same as in ?, but 4>* is a factor of 3 smaller (see Sect. l2.1.H l. 



Cluster Instrument and Reference Filter Exposure Time Depth (5tr, AB) Mask area 6' x 6', 2.2' x 2.2', 0.85' x 0.85' 

Abell 1689 ISAAC/VLT ? SZ 14.4 ksec 26l) 7% 12% 10% 

Abelll835 FORS/VLT ? V 19.6 ksec 27.0 6% 10% 8% 

AC114 SOFI/NTT? K 10.8 ksec 22.0 6% 9% 7% 

Table 2. Characteristics of the images used to produce the foreground object's mask for each cluster field, the three last columns 
representing the fractional surface area covered by the foreground galaxies in each cluster, for the three reference field of view, 
EMIR-like, JWST-like and NICMOS-like respectively 



they are centered on the critical lines in our study and avoid the 
crowded central regions of the cluster. In blank fields, this value 
doesn't exceed 3-4%. In the next sections, we have taken into ac- 
count this correction in the calculations of number counts, both 
in blank fields and lensing fields. 

Including the object mask M(x,y), number counts N(z,mo) 
are given by the following expression: 



N(z, m ) = cf> 



* J" M(x,y) 

x,y L(fi,z,mo) 

L(p, z, m ) 

■ exp " 



f Cv(x,y,z) ( L(p,z,m ) \ a 
J p(x,y,z) \ L* j 



^(z*) ^ x ^ 



(6) 



where p(x, y, z) is the magnification induced by the lensing 
field, C v {x,y,z) is the covolume associated with a single spatial 
resolution element pixel with Az = 1 and (L*, </>*, a ) are the 
parameters of the LF. 

A conservative upper limit of p — 25 was adopted in the 
vicinity of the critical lines in order to avoid singularities in 
the magnification/dilution determination. This is justified by the 
finite resolution of instruments, and the limited knowledge on 
the precise location of the critical lines at such high-z (typically 
~1"). 

When exploring the impact of cluster redshift, we assumed 
no evolution in the physical parameters r core , r cut and <x of in- 
dividual potentials, thus keeping the total mass of the cluster 
constant in this process. Variations in z c from the original red- 
shift ZcO °f the cluster produce a geometrical effect on the central 
positions (x ci , y ci ) of each PIEMD potential i, measured from a 
reference position (xo, yo) which coincides with the center of the 
BCG: 



%ci(.Zc) 

ya(zc) - yo 



Da(Zc) 

X = - - (X ci {Zc ) ~ X Q ) 



D A (z Co ) 

Da(Zc) 
D A (ZcS 



(ya(zc ) - yo) 



(7) 



Similarly, we produced fiducial masks at different cluster 
redshifts z c based on the reference mask at ZcO and adjust- 
ing the parameters of the galaxy ellipses, applying the same 
scaling relations on (x c ,y c ). The sizes a and b were scaled by 
D A (z c0 )/D A (z c ). 



3. Results 

As discussed above, lensing introduces two opposite trends on 
the observed sample as compared to blank fields : gravitational 
magnification by a factor p., increasing the expected number of 
sources and thus the total number of objects, and reduction of 
the effective surface by the same factor thus leading to a dilution 
in expected counts. This effect was first studied by ?. 

If we consider, for a given redshift z, the cumulative abun- 
dance of sources (per unit of solid angle) with a luminosity 
greater than L and by redshift bin, the magnification bias will 
change depending on p according to 



lensing 



(> L,z)=N(> L/p,z)/p(z) 
^/ (z) " 1 «(> L,z) 



(8) 
(9) 



where /3 is the logarithmic slope of n(L, z) assuming that this 
function is well represented by a power law in this interval of 
luminosities : f3 = -ii(lnn)/ii(lnL). The effect on number counts 
is as follows: 

- if P(z) > 1 the number counts will increase with respect to a 
blank field, and 

- if P(z) < 1 there will be an opposite trend: i.e. a depletion in 
number counts 

With increasing depth, the f3 parameter will decrease in a greater 
or lesser amount depending on the LF, the FOV (because it de- 
termines the mean p) and the redshift of sources, leading to a 
depletion in number counts in lensing fields as compared to a 
blank fields. With these simple considerations, we expect lens- 
ing clusters to be more efficient than blank fields in relatively 
shallow surveys. 
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The efficiency of using lensing clusters as gravitational tele- 
scopes to find high-z galaxies can be quantified with simple as- 
sumptions taking advantage of the properties of the sources ex- 
plained in Sect. [2] In this section, we discuss the results obtained 
by exploring the relevant intervals in the parameter space. We 
present a comparison between the number counts expected in 
lensing and blank fields, as a function of source redshift and for 
different LFs. The influence of lensing cluster properties and red- 
shift is also studied, as well as the expected field to field variance. 

3. 1 . The influence of the field of view 

Here we discuss on the influence of the FOV in the simulations 
for typical surveys. The influence of the limiting magnitude will 
be discussed in Sect. [4] Three different FOV are considered here: 

- 6' x 6' ("EMIR-like" aperture) 

- 2.2' x 2.2' ("JWST-like" or "WFC3/HST" aperture) 

- 0.85' x 0.85' (NICMOS/HST aperture) 

In the last case, the FOV is centered along the critical lines in 
order to achieve the highest mean magnification (Fig. [TJ. 

The limiting magnitude is AB <25.5, a value ranging be- 
tween L„(z=6) and 3L« - 5L* at redshift z s ~ 7 to 10. The cluster 
model corresponds to AO 14, but the results are qualitatively 
the same with other models. Fig. [3] displays the relative gain in 
number counts between lensing and blank fields as a function 
of sources redshift, for the three values of the FOV mentioned 
above, and for the three LF adopted in the present simulations. 

The largest gain is obtained for the smallest FOV, as expected 
from geometrical considerations, because the mean magnifica- 
tion decreases with increasing FOV, and in this case B > 1 given 
the shallow depth (A_B; 25.5). For a given FOV, the difference 
between lensing and blank field results strongly depends on the 
shape of the LF. Hence, the comparison between lensing and 
blank field number counts is likely to yield strong constraints on 
the LF, provided that field-to-field variance is small enough. This 
issue is addressed in Sect. 13.51 In the following subsections, we 
adopt a 6' x 6' FOV unless otherwise indicated. 

3.2. Lensing versus Blank field efficiency 

In this section, we study the effects of lensing clusters on source 
counts, using lensing models for the three reference clusters. 
We compute the expected number of sources brighter than mo, 
the typical apparent magnitude reached in ground-based near-IR 
surveys. The comparison between expected number counts of 
galaxies in a typical 6' x 6' FOV, up to mo < 25.5, per redshift 
bin Az - 1, in a blank field and in the field of a strong lensing 
cluster are presented in Fig.|4]in logarithmic scale. 

We also estimate the error on number counts due to the un- 
certainties on magnification factors (Sect. I2.1.2I I. The choice 
of the LF has no influence on the following results. Field to 
field variance dominate the error budget whatever the regime. 
Statistical errors and systematic errors on lensing models are 
smaller but not negligible as their contribution is less sensitive 
than field to field variance to the number of objects. In particu- 
lar, when the number of detected sources is relatively high (i.e. 
when field to field variance is relatively small), they reach ~ 15% 
of the error budget in the worst cases (e.g. for LF(a) at z ~ 6 in 
a 6' x 6' FOV), and typically < 2% when the number of sources 
is relatively small (e.g. for LF(c) at z ~ 8, for any FOV). 

As shown in Fig. |4] the presence of a strong lensing cluster 
has a dramatic effect on the observed number of sources, with 




LF(c) 



°6 8 10 12°6 8 10 12°6 8 10 12 

Z s 

Fig. 3. Relative gain in number counts between lensing and 
blank fields as a function of the source redshift, with Az - I, 
for different fields of view: 6' x 6' in black, 2.2' x 2.2' in blue 
dotted line and 0.85' x 0.85' in red dot-dashed line (from bottom 
to top). The three panels from left to right represent the 3 values 
of the LF, (a), (b) and (c) respectively. 



a positive magnification effect. Strong lensing fields are a factor 
between 2 and 10 more efficient than blank fields for the most 
optimistic LF (a), the gain increasing for the LFs (b) and (c), 
reaching a factor between 10 and 100 in the z ~ 6 - 12 domain. 
A positive magnification bias is observed, increasing with the 
redshift of the sources, and also increasing from optimistic to 
pessimistic values of the LF. This trend is indeed expected given 
the steep shape of the LF around the typical luminosity limits 
achieved in ground-based "shallow" surveys. 

Quantitatively (cf Tableland Fig.|4]i, if the LF for LBGs was 
nearly constant between z ~ 4 and 12, we could always detect 
at least one object over the redshift range of interest. At z ~ 6, 
we expect up to between 7-10 sources, and at z ~ 12 between 
0.7 and 1 galaxies should be detected in a lensing field. Even in 
a blank field, until z ~ 8 - 9 at least one LBG could be found in 
such a large field of view. With more realistic (pessimistic) val- 
ues of the LF (e.g. ??), blank fields are particularly inefficient as 
compared to lensing fields. The size of the surveyed area would 
need to increase by at least a factor of ~ 10 in order to reach 
a number of detections similar to the one achieved in a lensing 
field around z ~ 6 - 8, and this factor increases with redshift. 
Note however that given a limiting (apparent) magnitude, blank 
and lensing fields do not explore the same intrinsic luminosities 
(see also Sect©. 

As seen in Fig.|4]and Table[3] there are also some differences 
between the results found in the three lensing clusters, although 
they are smaller than the differences between lensing and blank 
fields for a given LF. The number of expected sources behind 
A 1689 is a factor of two (at z ~ 6) and a factor of three (z ~ 
8) larger than in the other clusters for the realistic LFs (b) and 
(c), whereas the difference is only ~ 10 - 30% for LF (a). The 
influence of lensing properties is studied in Sect. 13.41 

From the results above, it seems that lensing fields allow us 
to detect a larger number of z i 6 sources based on their UV con- 
tinuum, with some cluster to cluster differences. This result is 
essentially due to the shape of the LF. For magnitude limited 
samples selected within a given field of view, the positive mag- 
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nification bias increases with the redshift of the sources and de- 
creases with both the depth of the survey and the size of the sur- 
veyed area. The last trend is purely geometric, as discussed in the 
previous section. The differential behaviour between blank and 
lensing regimes strongly depends on the shape of the LF. The 
comparison between blank and lensing field observations could 
be of potential interest in constraining the LF, provided that field- 
to-field variance is sufficiently small. This issue is addressed in 
the following sections. 

3.3. Redshift of the lensing cluster 

The redshift of the lensing cluster is a potentially important pa- 
rameter when defining an "ideal" sample of gravitational tele- 
scopes. Based on geometrical considerations, we expect the 
magnification bias to decrease with cluster redshift (z c ) after 
reaching a maximum efficiency at some point, depending on 
cluster properties and the size of the surveyed field. The field 
of view considered here is typically a few square arcminutes, 
essentially including the region around the critical lines where 
magnification factors are the highest. Further down, we study 
the impact of z c on the magnification bias. 

Using the non-evolution assumption presented in Sect. 12.21 
we compute the expected number counts for the three refer- 
ence models (A1689, A1835 and AC1 14) with cluster redshifts 
ranging between z=0.1 and 0.8, with a Az = 0.1 step. A step 
Az = 0.05 is used in the z=0. 1-0.3 interval, in order to refine the 
sampling around the maximum value. We use the same depth 
and field size as in previous section. The effect of cluster red- 
shift is clearily seen in Fig. [^representing the number of objects 
as a function of cluster redshift (for the three reference models), 
at a fixed redshift of z=8 for sources. 

The global effect of z c on number counts as a function of the 
source redshift is displayed in Fig. [5] This figure directly com- 
pares to Fig. 2] in the previous section. Table [4] presents the z c 
value which corresponds to a maximum in the expected number 
counts at z=8. This value depends slightly on the source redshift 
and LF. In addition to the Az c — H — 0.05 when changing the 
LF, there is also an increase of z c with higher values of z s , up to 
+0.05 towards z s = 12. The search efficiency of distant galaxies 
in lensing fields is maximised when using clusters at low red- 
shift (z c ~ 0.1 - 0.3). Although the field of view considered here 
is relatively large for near-IR surveys and close to present-day 
cameras, it is the limiting factor at z c 5 0.1, where an increasing 
fraction of the strong-magnification area is lost with decreasing 
Zc- Also, in this z c regime, the field of view concentrates on the 
central region of the cluster where bright cluster galaxies mask 
an important fraction of the strong-magnification area. The high 
magnification region represents an increasingly small percent- 
age of the field with increasing z c - Number counts in this regime 
asymptotically tend towards a limiting value with increasing z c 
(Fig. 01, wich still is a substantial gain with respect to a blank 
field of the same size. The non-evolution assumption in cluster 
properties has a weak effect on this conclusion. Indeed, clusters 
far from relaxation will be even more ineffective as a gravita- 
tional lenses in the strongest magnification regime. The results 
obtained are hence an optimistic upper limit on number counts 
for realistic clusters beyond the optimal regime z c ~ 0.1 - 0.3. 

3.4. Influence of lensing cluster properties 

In this section we focus on the differences between lensing clus- 
ter properties and their influence on expected source counts. As 
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Fig. 6. Expected number of objects as a function of the cluster 
redshift for a fixed redshift of sources (z s - 8) and with the 
same depth and field size as in previous section. Three cluster 
models are displayed: A1689 (green dotted line), A1835 (blue 
dot-dashed line) and AC1 14 (red solid line). Panels from left to 
right display respectively LF (a), (b) and (c). 

Table 4. Redshift of the cluster which maximizes the number of 
objects detected at z=8 for the three LF respectively from top to 
bottom (a), (b) and (c) 



LF 


A 1689 


A1835 


AC114 


actual Zc 


0.184 


0.253 


0.310 


(a) 


0.10 


0.05 


0.25 


(b) 


0.15 


0.10 


0.25 


(c) 


0.20 


0.10 


0.30 



seen in previous sections, A1689-like clusters are expected to be 
more efficient irrespective of the cluster redshift. To understand 
this effect, we study the magnification regimes for a reference 
source plane fixed at z s = 8. The distribution of the magnifica- 
tion regimes in the image plane varies from cluster to cluster. 
Histograms in Fig.|7]represent the percentage of the image plane 
(for the 6' x 6' FOV) as a function of the associated magnifi- 
cation. To perform this calculation, cluster redshifts were stan- 
dardized to identical values for a better understanding of the phe- 
nomenon!. As seen in the figure, A 1689 shows a different regime 
at z c ~0.8 as compared to the other clusters. While the percent- 
age of the surface affected by strong magnification (fi > 10) does 
not exceed 5% in A1835 and AC1 14, it is as high as 8 - 10% in 
A1689, depending on z c . Nevertheless, this difference between 
clusters tends to fade with increasing cluster redshift due to pro- 
jection effects, the fraction of highly magnified pixels becoming 
smaller with respect to the whole FOV. We also note that AC1 14 
and A1835 models have a similar behaviour with minor differ- 
ences (A1835-like clusters being more effective at very low z c 
while the AC1 14 model is more efficient for intermediate z c ). 

Another way of understanding this phenomenom is pre- 
sented in the Fig. [8] where the effective covolume for the 6' X 6' 
FOV is traced as a function of the effective magnitude, for a 
magnitude-limited survey with AB < 25.5. Magnification in 
lensing fields provides an enhanced depth for a magnitude- 
limited survey, where the effective (lensing corrected) covolume 
surveyed decreases with increasing effective depth. The behav- 
ior of A1689 in Fig. [8] illustrates the situation for this particu- 
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Table 3. Total number of objects expected within a 6' x 6' FOV (up to Hab < 25.5, Az - 1) from the three LF adopted in these 
simulations. Uncertainties correspond to 1 - cr level in magnification and lensing modeling 
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Fig. 4. Comparison between the expected number counts of galaxies in a typical 6' x 6' FOV, up to Hab ^ 25.5, per redshift bin 
Az = 1, in a blank field (dashed lines) and in the field of lensing cluster (solid lines) (left to right respectively A1689, A1835 and 
AC1 14). Expected counts are obtained by the integration of 3 different luminosity functions (a), (b) and (c) from top to bottom. The 
limit of one source detected in the field of view is indicated by an horizontal line to guide the eye. 



larly efficient cluster, allowing us to study a ~ lOOMpc 3 vol- 
ume to AB ~ 29.0 with a relatively modest observational in- 
vestment. Except for some particularly efficient lensing clusters 
(such as A 1689), most lensing fields should behave the same 
way as A1835 or AC1 14. 

3.5. Field to field variance 

In this section, we address the expected field-to-field variance 
affecting our previous results in order to estimate its impact in 
blank and lensed fields. We used two different approaches: the 
two-point correlation function estimation proposed by ? and a 
pencil beam tracer through the Millenium simulation. 

The first estimate is based on the method implemented by ?. 
This method for the calculation of the cosmic variance is based 
on the two points correlation function £(r) of the sample (?). 
Field o field variance is given by 



f v f v d 3 x 1 d 3 x 2 {(\x 1 -X2\) 

J v j v d 3 ^d^ 2 



(10) 



where V represents the volume of the survey. 

We focus on the redshift interval z ~ 6 - 8 using the present 
"shallow" survey parameters (see Sect. 14. IK both in blank and 
lensing fields. Here we use the same parameters as in Sect. 13.21 
i.e. typical FOV 6' x 6', m A B = 25.5 and A z = 1. 
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Fig. 8. The effective (lensing corrected) covolume sampled at 
z=6.5-7.5 by each cluster is given as a function of effective Hab 
magnitude limit, for a magnitude-limited survey of 6' x 6' FOV 
with Hab ^ 25.5. The three clusters are displayed with the same 
colors and line codes as in Fig. [6] 



We define the total fractional error of the counts N following 
? (this is the so-called field-to-field standard deviation or the, 
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Fig. 5. The same as Fig. Hbut using different assumptions for the redshift of the clusters, with z c £ [0. 1 , 0.8] and a 0. 1 step (for more 
details see Sect. 13.31 ) 
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Fig. 7. Histogram representing the percentage of the surface (6' x 6' FOV) as a function of the magnification for different redshifts 
of cluster (z c =0.1, 0.2, 0.3, 0.8), using the same color codes for the three clusters as in Fig. |6](A1689: green dotted line, A1835: 
blue dot-dashed line and AC1 14: red solid line) 



again, improperly called "cosmic variance") as : 

V(jV 2 ) - (AQ 2 
v, = (11) 

<A0 

Results are presented in the Table [5]for the three LF consid- 
ered, and for the three typical clusters used in the simulations. 



We note an important field to field variance with such a lim- 
iting magnitude (AB ~ 25.5) either in blank or lensing fields 
due to the small number counts previously derived from calcula- 
tions (see Table [3j. Nevertheless, the variance is smaller behind 
gravitational telescopes with the same differential trends men- 
tioned before between the three clusters, i.e. A1689 exhibits a 
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stronger magnification bias (see Sect. 13. 4t than the other clus- 
ters which have a similar behavior. Besides, with increasing red- 
shift of sources, the expected number counts decrease leading to 
a larger field to field variance. 

The second estimate was based on the Millennium simula- 
tion, carried out by the Virgo Consortium and described in detail 
in ? and ?. The simulation follows N = 2160 3 particles of mass 
8.6 10 8 ft -1 M within a co-moving box of size 500 h 1 Mpc on a 
side. The cosmological model is a ACDM model with small dif- 
ferences in the cosmological parameters adopted in Sect. [1] but 
without impact on the final results. These cosmological parame- 
ters are consistent with recent determinations from the combined 
analysis of the 2dFGRS and 3rd year WMAP data (?). Given 
its high resolution and large volume, the Millennium simulation 
allows us to follow in enough details the formation history of 
a representative sample of high redshift galaxy environments. 
With these prescriptions and a realistic beam tracer we can study 
the field-to-field variations in the number counts of star forming 
galaxies at the epoch of interest. 

Our pencil beam tracer is similar to the one developed by ?. 
We trace through the simulation box a parallelepiped where the 
base is a parallelogram, whose size is given by the reference field 
of view in comoving units, and the depth is the comoving depth 
arbitrarily taken to Az = 1. The variation of angular distance 
versus redshift in the redshift interval of the selection window 
considered was properly taken into account. This 500/i~ l Mpc 
edge box is more than 2000 times larger than the effective vol- 
ume probed by the 6' x 6' FOV : ~ 10.7 x 10.7 x 323 .0(M pc / hf 
for instance at z=6. We carried out 10000 Monte Carlo realiza- 
tions of the beam-tracing procedure by randomly varying the ini- 
tial position of the beam in order to calculate the typical number 
counts of galaxies and the associated standard deviation in the 
field of view with the same hypothesis. 

Although this procedure is well suited to determine the field 
to field variance, several studies on this topic suggest an over- 
prediction on the abundance of massive galaxies at high redshift 
(e.g. ?). For this reason, we consider this second approach as a 
cross check yielding a lower limit for the field to field variance. 
Results obtained from the Millenium simulation are displayed in 
Table [6] They are in fair agreement with those obtained with the 
first method. 

Field to field variance on number counts obviously depends 
on the depth of the survey. In order to compare our results with 
existing photometric surveys, we calculated the number counts 
of sources in blank and lensing fields (here AO 14) with the 
evolving LF(c) for different deeper magnitude limits (m;, m 
27.0, 28.0 and 29.0), in our reference field of view using the 
same parameters as in Sect. 13. 2| (see Table |7j. The correlation 
function was used to derive the cosmic variance. The total frac- 
tional error (v r ) strongly decreases with increasing photometric 
depth, as expected given the increasing number of sources de- 
tected in such a large FOV (e.g. at m/,- m = 29.0, the total number 
of sources is 1000 times larger than in the "shallow" survey), 
both in blank and lensing fields. The fractional error appears 
slightly larger in lensing than in blank fields at z=6, but this 
effect reverses with increasing source redshift. These estimates 
for the blank field can be compared to present-day surveys. For 
instance, the field to field variations obtained by ? for a single 
ACS pointing at z ~ 6 for a limiting magnitude z%sq ~ 29 is 
35%. Using the same observational constraints (FOV, depth, ...), 
our simulations yield a v r ~ 30%, a value which is smaller but 
fairly compatible with the results quoted by ?. 



Table 6. Number counts for ot^b < 25.5 and field to field uncer- 
tainties (v f ) calculated from the Millenium simulation in a 6' x6' 
blank field, for different source redshifts 





z s = 6 


z s = 7 


Zs = 8 


(N) 


2.74 


1.10 


0.39 


V r 


68% 


100% 


168% 



4. Discussion 

4. 1 . Survey parameters and efficiency 

As discussed in Sect. I2.1.3I the FOV and the limiting magni- 
tude are two important survey parameters used in these simula- 
tions. The influence of the FOV for a fixed limiting magnitude 
strongly depends on the shape of the LF. The highest ratio in 
number counts between lensing and blank fields can be achieved 
with the smallest FOV due to simple geometrical considerations. 
This section specifically addresses the evolution on the survey 
efficiency in lensing and blank fields as a function of the limit- 
ing magnitude. 

For these purposes, we use the same approach as in Sect. l3.2l 
to derive number counts within a 6' x 6' FOV in blank fields 
and behind lensing clusters. AC1 14 is used here as a representa- 
tive lensing cluster. Fig.|9]displays the expected number counts 
as a function of the redshift of sources, for different depths 
(AB < 26.0,27.0,28.0 and 29.0). An opposite trend between 
blank and lensing fields appears, depending once again on the 
LF and on the redshift of sources. With increasing limiting mag- 
nitude, the efficiency of the survey towards a foreground clus- 
ter diminishes and becomes less effective than in blank fields 
leading to a negative magnification bias for the faintest limiting 
magnitudes (e.g. for LF(a) between AB ~ 27.0 - 28.0, for LF(b) 
between AB ~ 28.0 - 29.0 and for LF(c) beyond AB ~ 29.0). 
This trend, however, is highly sensitive to the FOV. In particular, 
the negative magnification bias appears towards the typical mag- 
nitudes achieved by space facilities (JWST). Fig. \W\ displays 
the same results as in Fig.|9]but for a 2.2' x 2.2' FOV (JWST- 
like). The main characteristics remain broadly unchanged, the 
general trends are just exacerbated, the inversion happening to 
lesser depth. 

Lensing and blank field surveys do not explore the same in- 
trinsic luminosities, as shown in Fig. QT| and Q~2] These figures 
compare the expected number density of sources as a function 
of their intrinsic UV luminosity (or equivalent SFR) for differ- 
ent limiting magnitudes ranging from AB;25.5 to 29.0. In the 
case of lensing fields, two different results are given, depending 
on the FOV around the cluster center. In this particular case, the 
source redshift is arbitrarily fixed to z s = 8, assuming a strongly 
evolving LF(c), and the lensing cluster is AC1 14. 

In summary, the number of z>8 sources expected at the typ- 
ical depth of JWST (AB ~ 28 - 29) is much higher in lensing 
than in blank fields if the UV LF is rapidly evolving with red- 
shift (LF(c)), as suggested by ?. The trend should be the oppo- 
site if the LF remains unchanged between z ~ 6 and 8. Lensing 
clusters are the only way to study the faintest building blocks 
of galaxies, with typical S FR <0.1 to lM /yr. On the contrary, 
wide field surveys covering 10 3 to 10 4 arcmirT 2 are needed to 
set reliable constraints on the brightest part of the LF at z J 6, i.e. 
for galaxies with SFR > l0M Q /yr. 
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Table 5. Number counts, field to field variance calculated with the correlation function both in blank and lensing fields, for z = 
6, 7 and 8 within a 6' x 6' FOV, for a shallow survey with AB < 25.5. Field to field variance in lensing fields includes 1 - cr 
magnification errors. 
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Table 7. Field to field variance for 3 different magnitude limits : m/,,„ = 27.0, 28.0 and 29.0, in a 6' x 6' blank field and lensing field 
(behind AC1 14) for the LF(c). Field to field variance in lensing fields includes 1 - cr magnification errors. 
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4.2. Influence of galaxy morphology and image sampling 



Gravitational magnification (e.g. in the tangential direction) 
induces an elongation of images along the shear direction 
while preserving the resolution in the perpendicular direc- 
tion and the surface brightness of high redshift galaxies. All 
the comparisons between lensing and blank fields in our sim- 
ulations assumed that observations were conducted with the 
same intrument setup in terms of FOV and spatial sampling, 
and with the same observational conditions, in particular the 
same limiting surface brightness and PSF. However, when 
comparing magnitude-limited samples in lensing and blank 
fields, it is worth discussing the influence of galaxy morphol- 
ogy and image sampling on the present results. In particu- 
lar, the evolution in the surface brightness of high redshift 
sources is susceptible to hinder the search efficiency in clus- 
ters if, for instance, number counts in clusters were domi- 
nated by sources below the limiting surface brightness. 



As explained in Sect. 12.1.11 all the previous results have 
been obtained assuming that galaxies at z > 7 are compact as 
compared to spatial sampling. Indeed, high redshift sources 
are expected to be very small, typically <0."10 on the sky, 
based on cosmological simulations (e.g. ?), in such a way that 
the high resolution capability of JWST is needed for resolv- 
ing such faint galaxies. Recent observations of LBGs candi- 
dates in the HUDF fully support this idea (???). In a recent 
paper, ? measured the average intrinsic size of z ~ 7-8 LBGs 
to be 0.7 + 0.3 kpc. These galaxies are found to be extremely 
compact, with very little evolution in their half-light radii 
between z ~ 6 and 7, roughly consistent with galaxies having 
constant comoving sizes, at least within the observed lumi- 
nosity domain ~ 0.1 - lL*(z = 3). Smaller physical sizes are 
expected for higher redshift and/or intrisically fainter galax- 
ies, based on the scaling of the dark matter halo mass or 
the disk circular velocity (?). This differential trend is actu- 
ally observed between the bright (~ 0.3 - 1 L*) and the faint 
(~ 0.12 - 0.3L*) samples of ?. 
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If all high-z galaxies exhibit the same compact and uni- 
form morphology, the effective mean surface brightness of 
a lensed galaxy will be brighter or fainter with respect to 
a blank field galaxy with the same apparent magnitude de- 
pending on the spatial resolution (in practice, the instrumen- 
tal PSF). The majority of lensed sources should remain spa- 
tially unresolved on their width on seeing-limited ground- 
based surveys, and even on their tangential direction up to 
a gravitational magnification \i ~ 5 - 10. Hence, the appar- 
ent surface brightness of a lensed source is actually brighter 
than that of a blank field galaxy of similar apparent magni- 
tude (by roughly -2.5log(p.) mags for a spatially unresolved 
galaxy). This situation is typically found in the "shallow and 
wide" near-IR surveys discussed above (e.g. for the 6' x 6' 
FOV), where lensing clusters are particularly efficient. 

On the contrary, for a fixed apparent magnitude, the 
effective mean surface brightness of a lensed galaxy is ex- 
pected to become fainter with respect to a blank field galaxy 
when the image resolution is similar or better than its (lensed 
maximum) half-light radius, reaching ~ 2.5log(/j.) mags in 
the worst case. This situation is typically expected in the 
"deep and narrow" near-IR surveys with space facilities. 
In practice, the best spatial resolution presently achieved 
with HST/WFC3 in the near-IR is ~ O. 'IO, reaching 0."065 
with JWST/NIRCam, i.e. the typical size of the brightest 
z ~ 7 - 8 LBGs candidates presently identified. Therefore, 
the majority of lensed sources should remain spatially unre- 
solved on their width. A lensed source entering the apparent- 
magnitude limited sample because of its magnification yu has 
also a smaller physical size, by a factor of /j (assuming a con- 
stant M/L scaling with the halo mass) or /j 1/2 (assuming a 
constant M/L scaling with the halo circular velocity), leading 
to an apparent increase on its surface brightness with respect 
to blank-field observations of the same galaxy. Given the spa- 
tial resolution achieved with HST and JWST, this intrinsic- 
size effect tends to compensate the image dilution described 
above, in such a way that the actual surface brightness of the 
lensed galaxy should get close to the surface brightness of a 
blank field galaxy of similar apparent magnitude. 

For the reasons explained above, and to the best 
of present knowledge, we do not expect the apparent- 
magnitude limited number counts derived in clusters to be 
strongly biased by sources below the limiting surface bright- 
ness, provided that the usual scalings apply to the size of 
high-z sources. 

4.3. Comparison with current survey results 

We have compared our simulation results to recent observa- 
tions looking for high-z LBGs. For instance, the discovery of 
a bright z = 7.6 ± 0.4 lensed galaxy by ?, with AB=24.7 (intrin- 
sic AB ~ -22.6), in a 2.5' x 2.5' FOV survey around A1689 is 
in fair agreement with our expectations. Indeed, given the survey 
characteristics and including ~ 100% variance for z s = 7.6 + 0.4, 
we expect between 0.2 and 0.8 such bright objects in this lensing 
field, if the LF remains constant between z ~ 4 and 8 (LF(a)). 
In case of a strongly evolving LF(c), the expected number of 
sources in this survey is 0. 12 (i.e. ranging between and 0.5 with 
~ 200% variance) making the discovery of this bright source 
particularly fortunate. Our results for lensing fields are also con- 
sistent with the number of z ~ 7.5 LBGs found by ?, to the depth 
of their survey, using LF(b) or (c). Quantitatively, ? detected 5 
sources with 12 pointings over 6 clusters. With our simulations, 
2.6*9 -I objects with a variance of ~ 75% are expected with the 




Fig. 9. Expected number counts of objects as a function of the 
redshift of sources in a 6' x 6' FOV, for different limiting magni- 
tudes 26.0, 27.0, 28.0 and 29.0 from bottom to top respectively. 
This calculation is provided both in blank (dotted line) and lens- 
ing fields (solid line) (here AC1 14) and for the three LFs (from 
right to left, (a) in red, (b) in blue and (c) in black) 




Fig. 10. The same as Fig.|9]but for a 2.2' x 2.2' FOV. Some dif- 
ferences appear in comparison with the Fig. [9] For example, the 
total numbers of high-z galaxies expected behind lensing clus- 
ters (solid lines) and the field (dashed lines) are much larger at 
low limiting magnitude (juab = 26.0) but this phenomenon! is 
reversed for deeper surveys {mab = 28.0 - 29.0) (see text for 
details). 



LF(c) model. We also compared with the surface density of z ~ 7 
candidates in the deep near-IR data behind clusters obtained by 
?. ? found a surface density of 0.05*911 arcmirC 1 with AB < 25.5 
with a typical NICMOS3 FOV. With the strongly evolving LF(c) 
and the same survey characteristics used in our simulations, we 
expect a surface density of 0.01 arcmirC 1 behind a typical clus- 
ter such as AC1 14, with a variance of ~ 1 10%. This result shows 
a relatively good agreement by taking into account field to field 
variance. 
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Fig. 11. Cumulative surface density of sources as a function of 
their intrinsic UV luminosity, in blank fields (blue solid line) 
and in the lensing fields with FOV ~ 5 arcmin 2 (JWST-like, red 
solid line), for different photometric depths ranging from shal- 
low (AB ~ 25.5) to deep (AB ~ 29.0) surveys and a strongly 
evolving LF(c). The source redshift is arbitrarily fixed at z=8, 
with Az = 1 
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Fig. 12. Cumulative surface density of sources as a function of 
their intrinsic UV luminosity and SFR, in blank fields (blue solid 
line) and in the lensing fields with FOV 6' x 6' (red dotted 
line) and 2.2' x 2.2' (JWST-like, red dashed line), for differ- 
ent photometric depths ranging from shallow (AB ~ 25.5) to 
deep (AB ~ 29.0) surveys and a strongly evolving LF(c). The 
mean magnification over the whole the field is used to derive 
the lensing points, the true distribution is displayed in Fig. QT| 
The source redshift is arbitrarily fixed at z=8, with Az = 1. 
The conversion from absolute magnitude to SFR is provided in 
Sect. l2.1.3l using the calibrations from ?. 



4.4. Lyman Break versus NB searches 

In this section, we discuss on the relative efficiency of blank and 
lensing fields on the detection of LAEs based on either NB sur- 
veys or the spectroscopic follow up of LBGs at z>6. Although 
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Fig. 13. Cumulative surface density of observed sources as a 
function of their Lyman alpha luminosity (6' x 6' FOV, A z = 1, 
redshift of sources fixed at 6.6 for the LF(c)). The density is cal- 
culated in blank (blue solid line) and in lensing fields (red solid 
line) for different limiting magnitudes (from right to left: 25.5, 
26.0, 27.0, 28.0 and 29.0). Dashed lines display number counts 
corrected by transmission value e ~ 10% (dashed red line in 
lensing fields and dashed blue line in blank field). For compari- 
son, raw number counts extracted from the spectroscopic sample 
of LAEs by ? are also given (black solid line). The number den- 
sity derived from ? at z = 6.96 is also indicated, together with 
corresponding error bars (see text for details). As in Fig.[l2j the 
magnification used to derive the lensing points is averaged over 
the entire field. 



the observational effort required to select z > 7 candidates using 
the dropout technique seems relatively cheap as compared to the 
NB approach, the two approches are complementary, as empha- 
sized by the fact that many objects found by Lyman a emission 
remain weak or undetected in the continuum (e.g. ?, ?, ?, ?). 
A quantitative comparison between the properties of LAEs and 
LBGs at z > 7 within the same volume should provide impor- 
tant information on the Lyman a transmission, SFR and other 
properties of these high-z galaxies. 

Since the pioneering Large Area Lyman Alpha Survey 
(LALA, ?, ?), different NB surveys in blank fields have pro- 
vided interesting galaxy samples in the z~ 5 - 7 interval, e.g. 
the large sample of Lya emitters at z ~ 5.7 by ?, the z=6. 17 and 
6.53 galaxies found respectively by ? and ?, the two z ~ 6.6 
galaxies detected by ?, and the galaxy at a redshift z=6.96 
found by ?. In the latter case, which should be representative 
of z~ 7 samples, the authors used a combination of NB imaging 
at 8150A (SuprimeCam) and broad-band photometry in the op- 
tical bands to select candidates for a subsequent spectroscopic 
follow up with DEIMOS/Keck. Their confirmation rate is rel- 
atively high (18 sources out of 26 candidates), leading to 0.03 
sources/arcmin 2 and redshift bin 5z = 0.1. Similar results are 
reported by ?. All these sources have important Lyman a fluxes 
(a few xl0~ 17 erg cirT 2 s~'), and display broad Lyman a lines 
(cr,, ~ 200 km/s). A strong evolution is found in the number den- 
sity of LAEs at z > 7 with respect to the z~ 5 - 7 interval (?, ?, 
?)■ 
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The number of LAEs expected within a sample of LBGs at 
z ~ 6 can be estimated using the distribution of Lyman-ff equiv- 
alent widths derived for the spectroscopic sample of LBGs at 
Z ~ 3 by ?, assuming no evolution in the population of LAEs 
with respect to LBGs. This simplistic scaling should be enough 
for the simulation needs. We introduce a factor e, defined below, 
which can be linked to the Lyman a transmission as follows: 



where L1500 is the UV monochromatic luminosity at 1500 A, 
LLya is the Lyman-ff luminosity and WLya is the Lyman-c equiv- 
alent width. With this simple assumption, the average value 
for the Lyman-a equivalent width is ~ 10 A, corresponding to 
e ~ 10%. This value can be used to derive a rough estimate of 
expected number density of LAEs, from a population of LBGs 
In addition, the number density is also corrected to take into 
account of the fraction of the LBGs sample displaying Lya in 
emission. 

Fig. Q~3] displays the cumulative number counts of sources 
at z ~ 6.6 integrated from the LF(c) as a function of the 
Lyman-a luminosity, scaled according to the UV luminosities 
(cf Sect. I2.1.U in the typical 6' x 6' FOV, together with a com- 
parison of observations in a similar redshift domain (z s ~ 6.6 for 
Kashikawa spectroscopic sample of LAEs and ? at z=6.96). 

The number density of LBGs at z=7 (with 8z = 0.1, 
close to the band-width of NB surveys) ranges between 0.001 
(LF(c)) and 0.02 (LF(a)) sources/arcm/n 2 for a survey limited 
to H(AB) <25.5, depending on the LF. Lensing clusters improve 
these numbers by a factor ranging between 6 (for LF(c)) and 2 
(for LF(a)). In case of a deep survey limited to H(AB) < 29.0, the 
number densities reach 1 (LF(c)) to 2 (LF(a)) sources/ 'a rcmin 2 . 
In this case, there is a negative magnification bias of the order 
of 20%. These numbers, obtained with a simplistic model, are 
between a factor of ~ 10 (for bright sources) and a few (for 
faint sources) smaller than the number densities obtained by ? 
for their spectroscopic sample. With increasing z s (see Fig. |9j 
for instance at z=9 with the strongly evolving LF(c), no sources 
can be detected for a shallow survey limited to H(AB)<25.5 
and for a deeper limited survey (H{AB) <29.0), a minimum of 
3 arcmin 2 surveyed area is needed to obtain 1 source in a blank 
field. In a lensing field with the (LF(c)), these number densi- 
ties reach 0.002 for H(AB) <25 .5 and 0.32 sources/arcm/n 2 for 
H(AB) $29.0. The relatively low-efficiency of lensing clusters 
with NB techniques in the z > 9 domain has been recently con- 
firmed by the results of ?. 

The preselection of z > 6 - 7 candidates in lensing fields has 
two main advantages with respect to blank fields. In the shallow 
(AB < 25.5) regime, there is an increase by a factor ~ 8 - 10 on 
the number of sources detected and a moderate gain in depth for 
a given exposure time (i.e. ~ 0.5 magnitudes at AB ~ 25.5). In 
the deep-survey regime {AB < 28 -29), there is a gain in intrinsic 
depth, for a number of candidates which remains essentially con- 
stant (i.e. ~ 0.8 mags gain at AB < 28). The relative efficiency 
of lensing with respect to blank field counts in Fig. [T2l depends 
on the FOV. The two predictions get close to each other with 
increasing values of the FOV in lensing surveys, and the trend 
goes in the opposite direction for smaller FOV. This trend is the 
same for both LBGs and LAEs. To explore the bright end of the 
LF, blank field surveys are needed with a large FOV, whereas 
lensing clusters are particularly useful to explore the faint end of 
the LF. This trend is further discussed in the next Section. 



4.5. Towards the ideal survey: constraining the Luminosity 
Function of high-z sources 

All present photometric surveys aimed at constraining the UV 
LF at z>7, either space or ground-based, are still dramatically 
small in terms of effective surface. Wide and deep optical+ 
near-IR surveys in lensing and blank fields are needed to set 
strong constraints on the LF and on the star-formation density 
at z>7. An important issue is the combination between photo- 
metric depth and surveyed area which is needed to identify a 
representative number of photometric candidates, or to reach a 
significant non-detection limit in order to constrain the LF of 
z>7 sources. 

There are three different aspects to consider when design- 
ing an "ideal" survey aiming at constraining the LF: the depth 
and the area of the survey, and the corresponding field to field 
variance. In order to address these issues, we have computed 
the expected field to field variance corresponding to lensing and 
blank field surveys, for different survey configurations (area and 
depth). A summary of these results is given in Table[8]for differ- 
ent number of lensing clusters, and for two representative depths 
in the H-band (i.e. a "shallow" survey with AB < 25.5, and a 
"deep" survey with AB < 28.0) assuming a strongly evolving 
LF(c) in all cases. This table complements the results given in 
Tables|5]and|7]for blank and lensing fields as a function of depth. 
In all cases, we use AC1 14 as a reference for lensing clusters. 

Regarding field-to-field variance in number counts, results 
are expected to be similar in blank and lensing fields for a rela- 
tively wide FOV (~ 40 - 50 arcmin 2 ; see Sect. l3.5l and Tabled. 
As shown in Table [8] a deep lensing survey using ~ 10 clusters 
should be able to reach a variance < 20% on sources at 6 < z < 8, 
irrespective of the actual LF. This value is better than present-day 
photometric surveys in blank fields, typically reaching 30-35% 
for AB < 29.0 (e.g. ?), which in turn is rather close to what could 
be achieved in a single lensing cluster for AB < 28.0. 

A different survey strategy consists of increasing the number 
of lensing clusters with a shallow limiting magnitude. In this 
case, a few tens of lensing clusters (typically between 10 and 
50, depending on the LF) are needed to reach a variance of ~ 
30% at z< 8. Note that the difference in exposure time between 
the shallow and deep surveys reported in Table [8] is a factor ~ 
100, and that ~ 10 pointings are needed on blank fields in order 
to reach the same number of z ~ 6 - 8 sources as in a single 
"shallow" lensing field. 

In the case of a strongly evolving LF(c), photometric surveys 
should reach a minimum depth of AB ~ 28 to achieve fair statis- 
tics on z ~ 9-10 sources using a lensing cluster (Table|7]i. In this 
case we expect between 20 (z=9, ~ 40% variance) and 8 (z=10, 
~ 40% variance) sources per lensing cluster in a ~ 40 arcmin 2 
FOV. The efficiency is a factor of 10 smaller at z ~ 10 in blank 
fields. Fair statistics at z ~ 12 should require a minimum depth 
close to AB ~ 29 both in lensing and blank fields. 

Constraining the LF of star forming galaxies at z > 7 should 
require the combination of blank and lensing field observations. 
This is illustrated for instance in Fig.fTTIand [T2]for an example 
at z s — 8. A survey towards a lensing cluster has several advan- 
tages. It increases the total number of sources available for spec- 
troscopic follow up, and it helps extending the sample towards 
the faint edge of the LF and towards the highest possible limits 
in redshift. On the other hand, blank fields are needed to achieve 
fair statistics on the bright edge of the LF. Thus an "ideal" survey 
should combine both blank and lensing fields. Given the num- 
bers presented in previous sections, a blank field used for these 
purposes should be a factor ranging between ~ 10 and 100 times 
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larger than a lensing field (depending on the redshift domain, 
photometric depth, and actual LF) in order to efficiently com- 
plete the survey towards L > L* . This should be possible with 
the new upcoming surveys, such as the WIRCAM ultra deep 
survey (WUDS) at CFHT (~ 400 arcmin 2 , with YJHK <25.5), 
UKIDSS-UDS (~ 2700 arcmin 2 , with YJHK< 25) or Ultra- 
Vista (~ 2600 arcmin 2 , with YJH < 26, K< 25.6). The optimum 
number of lensing fields ranges between ~ 10-20 (for z ~ 6 - 8 
studies with "shallow" photometry) to a few (for "deep" surveys 
targeting z ~ 8 - 12 sources). 

5. Summary and conclusions 

We have evaluated the relative efficiency of lensing clusters with 
respect to blank fields in the identification and study of z > 6 
galaxies. The main conclusions of this study are given below. 

For magnitude-limited samples of LBGs at z > 6, the magni- 
fication bias increases with the redshift of sources and decreases 
with both the depth of the survey and the size of the surveyed 
area. Given the typical near-IR FOV in lensing fields, the max- 
imum efficiency is reached for clusters at z ~ 0.1 - 0.3, with 
maximum cluster-to-cluster differences ranging between 30 and 
50% in number counts, depending on the redshift of sources and 
the LF. 

The relative efficiency of lensing with respect to blank fields 
strongly depends on the shape of the LF, for a given photometric 
depth and FOV. The comparison between lensing and blank field 
number counts is likely to yield strong constraints on the LF. 

The presence of a strong-lensing cluster along the line of 
sight has a dramatic effect on the observed number of sources, 
with a positive magnification effect in typical ground-based 
"shallow" surveys (AB < 25.5). The postive magnification bias 
increases with the redshift of sources, and also from optimistic 
to pessimistic values of the LF. In case of a strongly evolving 
LF at z > 7, as proposed by ?, blank fields are particularly in- 
efficient as compared to lensing fields. For instance, the size of 
the surveyed area in ground-based observations would need to 
increase by a factor of ~ 10 in blank fields with respect to a 
typical ~ 30-40 arcmin 2 survey in a lensing field, in order to 
reach the same number of detections at z ~ 6 - 8, and this merit 
factor increases with redshift. All these results have been ob- 
tained assuming that number counts derived in clusters are 
not dominated by sources below the limiting surface bright- 
ness of observations, which in turn depends on the reliability 
of the usual scalings applied to the size of high-z sources. 

Ground-based "shallow" surveys are dominated by field-to- 
field variance reaching ~ 30 to 50% in number counts between 
z~6 and 8 in a unique ~ 30 - 40 arcmin 2 lensing field survey (or 
in a 400 arcmin 2 blank field), assuming a strongly evolving LF. 

The number of z>8 sources expected at the typical depth of 
JWST (AB ~ 28 - 29) is much higher in lensing than in blank 
fields if the UV LF is rapidly evolving with redshift (i.e. a factor 
of ~ lOatz ~ 10 withAfi<28). 

Blank field surveys with a large FOV are needed to probe the 
bright edge of the LF at z > 6 - 7, whereas lensing clusters are 
particularly useful to explore the mid to faint end of the LF. 
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Table 8. Field to field variance, including 1 — cr errors on the magnification factor, expected in a lensing survey, as a function 
of the number of clusters, for the three LFs (a), (b) and (c) (from top to bottom respectively). Simulations were performed using 
the same approaches in Sect. [33] (i.e. 6' x 6' FOV, ...) for two different depths: AB < 25.5 (shallow survey) and AB < 28.0 (deep 
survey) 
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